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Abstract Active Galactic Nuclei (AGN) represent the growth phases of the super- 
massive black holes in the center of almost every galaxy. Powerful, highly ionized 
winds, with velocities ~ 0.1 — 0.2c are a common feature in X-ray spectra of lu¬ 
minous AGN, offering a plausible physical origin for the well known connections 
between the hole and properties of its host. Observability constraints suggest that 
the winds must be episodic, and detectable only for a few percent of their lifetimes. 
The most powerful wind feedback, establishing the M — a relation, is probably not 
directly observable at all. The M — a relation signals a global change in the nature 
of AGN feedback. At black hole masses below M — a feedback is confined to the 
immediate vicinity of the hole. At the M — a mass it becomes much more energetic 
and widespread, and can drive away much of the bulge gas as a fast molecular 
outflow. 

Keywords Supermassive black holes, accretion, M — a relation, X-ray winds, 
molecular outflows, quenching of star formation 


1 INTRODUCTION 

1.1 SMBH Scaling relations 

Astronomers now generally agree that the center of almost every galaxy but the 
smallest contains a supermassive black hole (SMBH). In recent years it has become 
clear that the mass M of the hole correlates strongly with physical properties of the 
host galaxy. In particular the hole mass M appears always to be a fairly constant 
fraction of the stellar bulge mass M;,, i.e. 

M ~ (1) 

(Haring & Rix, 2004). Even more remarkably, observations give a tight relation of 
the form 

M ~ 3 X lO®M0cr^oo (2) 
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between the SMBH mass and the velocity dispersion a = 200cr2oo kms“^ of the 
host galaxy’s central bulge, with a ~ 4.4±0.3 (Ferrarese & Merritt 2000, Gebhardt 
et al. 2000; see Kormendy & Ho 2013 for a recent review). For many practical cases 
the relation is more conveniently written as 

M ~ 2 X lo'^M0crfoo (3) 


with now cr = lOOaioo kms~^. 

Since observationally determining the SMBH mass generally involves resolving 
its sphere of influence, of radius 
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with M = IO^MsMq = lO^MyM©, Q may represent a maximum SMBH mass for 
a given velocity dispersion a (Batcheldor 2010). 


1.2 Binding Energies 

At first sight the relations (HU appear surprising. For dl shows that the black 
hole’s gravity has a completely negligible effect on its host galaxy, which in most 
ways must be quite unaware of its existence. But we know (Soltan 1982) that the 
hole grew largely as a result of luminous accretion of gas. This released energy 

Ebh ~ 2 X 10® Vg erg, (5) 

where r] 0.1 is the accretion efficiency, far larger than the binding energy 

-^bulge ~ 8 X lO^^Mgaioo erg (6) 

of a host bulge of stellar mass ~ 10®M. 

The vast difference in these two numbers suggests that the host must notice 
the presence of the hole through its energy output, even though it is utterly in¬ 
significant in all other ways. We can already see how the black hole mass might 
correlate with galaxy properties - the hole grows by accreting gas, but in doing 
this communicates some of its huge binding energy Ebh back to the gas reservoir, 
and so potentially limits its own growth. This suggests that the most relevant 
quantity to compare with Ebh is not Ebuigei but instead the gravitational binding 
energy of the bulge gas alone, i.e. 

F/gas = fg Ehulge (7) 

where fg<lis the gas fraction. In the following we take this as fg ~ 0.16, the 
cosmological mean value, giving a typical relation 

F-bh ~ 2000Figas (8) 

for a black hole close to the M — a relation (the rhs has an implicit factor ~ 
^200/^s)- 

This picture requires the black hole to communicate some of its accretion 
energy to its host. But this process cannot be very efficient, as otherwise the hole 
could disrupt the host entirely, or at the very least remove a large fraction of its 
gas. In this sense, the galaxy bulge leads a precarious existence. For much of its 
life it can ignore the threat that the SMBH poses, but we will see that in the end 
this is always decisive if accretion continues. 
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1.3 Communicating the Energy: Feedback 

There are two main ways that the SMBH binding energy can potentially inter¬ 
act with its surroundings. By far the larger is direct radiation: after all, this is 
how all the accretion energy is initially released. But we know from observation 
that most light escapes relatively freely from active galactic nuclei (AGN). This 
suggests that radiation is in general not the main way the SMBH affects its host, 
and we will discuss in detail why this is so in Section 7.4. The second form of 
coupling SMBH binding energy to a host bulge is mechanical. The huge SMBH 
accretion luminosity drives powerful gas flows into the host, making collisions and 
communication inevitable. One form of flow often mentioned is jets — highly col¬ 
limated flows driven from the immediate vicinity of the SMBH (see Fabian 2012 
for a review). To turn these into a way of affecting most of the bulge requires a 
way of making the interaction relatively isotropic, perhaps with changes of the jet 
direction over time. Here we will mainly consider another form of mechanical com¬ 
munication which automatically has this property already. This is the observed 
presence in many AGN of near-isotropic winds carrying large momentum fluxes. 


1.4 Powerful ionized winds 

Early X-ray observations of AGN yielded soft X-ray spectra frequently showing 
the imprint of absorption from ionized gas, the ‘warm absorber’; hereafter WA 
(Halpern 1984, Reynolds & Fabian 1995). More recent observations have found 
at least 50% of radio-quiet AGN showing WAs in their soft X-ray (~0.3-2 keV) 
spectra. The limited spectral resolution of the Einstein Observatory and ASCA 
observations prevented important parameters of the WAs, in particular the out¬ 
flow velocity and mass rate, to be determined with useful precision. The higher 
resolution and high throughput afforded by contemporary X-ray observatories, 
Chandra, XMM-Newton and Suzaku has transformed that situation over the past 
decade, with the WA being shown, typically, to be dominated by K-shell ions of the 
lighter metals (C, N, O, etc) and Fe-L, with outflow velocities of several hundred 
km s“^ (Blustin et al. 2005, McKernan et al. 2007). 

A more dramatic discovery made possible with the new observing capabilities 
was the detection of blue-shifted X-ray absorption lines in the iron K band, in¬ 
dicating the presence of highly ionized outflows with velocities v ~ 0.1 — 0.25c 
(Ghartas et al. 2002; Pounds et al. 2003; Reeves et al. 2003). In addition to adding 
an important dimension to AGN accretion studies, the mechanical power of such 
winds, which for a radial flow depends on v^, was quickly recognized to have a 
wider potential importance in galaxy feedback. 

Additional detections of high velocity AGN winds were delayed by the low 
absorption cross section of such highly ionized gas, combined with strongly blue- 
shifted lines in low-redshift objects coinciding with falling telescope sensitivity 
above ~ 7 keV. However, further extended observations, particularly with XMM- 
Newton, found evidence in 5 additional AGN for outflow velocities of ~ 0.1 — 
0.2c (Gappi et al. 2006). Some doubts remained as to how common high velocity 
outflows were, as the majority of detections were of a single absorption line (with 
consequent uncertainty of identihcation), and had moderate statistical significance, 
raising concerns of ‘publication bias’ (Vaughan and Uttley 2008). In addition, only 
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for PG1211+143 had a wide angle outflow been directly measured, confirming a 
high mass-rate and mechanical energy in that case (Pounds & Reeves 2007, 2009) 

These residual doubts were finally removed following a blind search of extended 
AGN observations in the XMM-Newton archive (Tombesi et al. 2010), finding com¬ 
pelling evidence in 13 (of 42) radio quiet objects for blue-shifted iron K absorption 
lines, with implied outflow velocities of ~ 0.03 — 0.3c. A later search of the Suzaku 
data archive yielded a further group of strong detections, with a median outflow 
velocity again ~ 0.1c (Gofford et 2013). In addition to confirming that high veloc¬ 
ity, highly ionized AGN winds are common, the yield from these archival searches 
shows the flows must typically have a large covering factor, and therefore be likely 
to involve substantial mass and energy fluxes. 

The observed distributions of velocity, ionization parameter and column den¬ 
sity are compatible with Eddington winds launched from close to the black hole, 
where the optical depth res ~ 1, and carrying the local escape velocity (King & 
Pounds 2003). However, as the mean luminosity in most low-redshift AGN is on 
average sub-Eddington, such winds are likely to be intermittent, a view supported 
by repeated observations and by the range of observed column densities. 

For the best-quantified high-velocity outflow (the luminous Seyfert PG1211-1-143), 
in which a wide-angle flow was directly measured (Pounds & Reeves 2007, 2009), 
the wind appeared to have more energy than needed to unbind the likely gas mass 
of the observed stellar bulge. This suggested that the energy coupling of wind to 
bulge gas must be inefficient, as seen in the discussion following equation (|8l). Ev¬ 
idence that the fast wind in NGC 4051 is shocked at a distance of ~ 0.1 pc from 
the black hole offers an explanation of why such powerful winds do not disrupt 
the bulge gas: strong Compton cooling by the AGN radiation field removes most 
of the wind energy before it can be communicated. 


2 THE OBSERVATIONAL EVIDENCE FOR ULTRA FAST OUTFLOWS 

As noted above, the requirement of X-ray observations with high sensitivity and 
good spectral resolution over a wide energy band delayed the discovery of powerful, 
highly ionized winds from non-BAL AGN until the launch of Chandra and XMM- 
Newton. A decade after the first reports (Pounds et al. 2003, Reeves et al. 2003), 
high-velocity (v~0.1c), highly-ionized winds are now established to be common in 
low redshift AGN. 


2.1 The fast outflow in PG1211-I-143 

Exploring the nature of the ’soft excess’ in a sample of luminous Palomar Green 
AGN was a primary target in the Guaranteed Time programme awarded to Martin 
Turner, Project Scientist for the EPIC Camera on XMM-Newton (Turner et al. 
2001). At that time the X-ray spectrum in AGN above ~1 keV was expected to 
be a rather featureless power law apart from a fluorescent emission line at ~6.4 
keV from near-neutral Fe. One source, PG1211-I-143, showed a surprisingly ’noisy’ 
X-ray spectrum which one of us (KP) volunteered to explore. 

PG1211-I-143, at a redshift of 0.0809 (Marziani et al. 1996), is one of the 
brightest AGN at soft X-ray energies. It was classified (Kaspi et al 2000) as a 



AGN Outlflows and Feedback 


5 


Narrow Line Seyfert 1 galaxy (FWHM H/3 1800 km s“^), with black hole mass ~ 
4 X 1O^M0 and bolometric luminosity 4 x 10^®erg s“^, indicating a mean accretion 
rate close to Eddington. 

Analysis of the unusual spectral structure in the 2001 XMM-Newton observa¬ 
tion of PG1211-I-143 showed it to be dominated by blue-shifted absorption lines 
of highly ionized metals, providing the first evidence for a high velocity ionized 
outflow in a non-BAL AGN, with the initial identication of a deep blue-shifted Fe 
Lyman-a absorption line indicating an outflow velocity of ~ 0.09c (Pounds et al. 
2003). That observation, closely followed by the detection of a still higher outflow 
velocity from the luminous QSO PDS 456 (Reeves et al. 2003), attracted wide at¬ 
tention, potentially involving the ejection of a significant fraction of the bolometric 
luminosity, and perhaps characteristic of AGN accreting near the Eddington rate 
(King & Pounds 2003). 

Appropriately for such an unexpected discovery, the validity of the high veloc¬ 
ity in PG1211-I-143 was not unchallenged. The near-coincidence of the observed 
absorption line blueshift and the redshift of the host galaxy was a concern, notwith¬ 
standing the uncomfortably high column density of heavy metals implied by a local 
origin. Then, in a detailed modelling of the soft X-ray RGS data, Kaspi & Behar 
(2006) found only a much lower velocity. Any doubts relating to the absorption 
being local were removed, however, by a revised velocity of 0.13 — 0.15c based on 
the inclusion of additional absorption lines from intermediate-mass ions (Pounds 
& Page 2006), and when repeated observations of PG1211-I-143 demonstrated that 
the strong Ee K absorption line was variable over several years (Reeves et al. 2008). 

Here we use the 2001 XMM-Newton observation of PG1211-I-143 with the pn 
camera (Strueder et al. 2001) to illustrate the two methods used then - and since 
- to parameterise the ionized outflow. Figure 1 shows the ratio of EPIG pn data 
to a simple power law continuum, with a deep absorption line seen near 7 keV and 
additional spectral structure at ~ 1 — 4 keV. 

Fitting a negative Gaussian to the deep ~ 7keV absorption line (figure 2, top 
panel), finds an observed line energy of 7.06 ± 0.02 keV, or 7.63 ± 0.02keV at the 
AGN redshift of 0.0809. The line is clearly resolved, with la width ~ 100 ± 30 
eV. Assuming identification with the Fe XXV resonance (6.70 keV rest energy), 
the blueshifted line corresponds to an AGN outflow velocity v~0.122±0.005c. The 
most likely a priori alternative identification, with the Fe XXVI Lyman-a line 
(6.97 keV rest energy), conservatively adopted in the initial analysis (Pounds et 
al. 2003), yields a lower outflow velocity v ~ 0.095 ± 0.005c. 

An alternative procedure, which also provides additional parameters of the gas 
flow, requires full spectral modelling, as in Pounds & Page (2006), and more widely 
in recent outflow studies (Section 2.3). For the 2001 XMM-Newton pn spectrum 
of PG1211-I-143, modelling the absorption from 1-10 keV with a photoionized gas 
derived from the XSTAR code of Kallman et al (1996) gives an excellent fit, for 
a column density N^ ~ 3.2 ± 0.7 x 10^^ cm“^, ionization parameter log^ = 2.7 ± 
O.lergcms”^, and an outflow velocity (in the AGN rest frame) v ~ 0.149±0.003c. 
The model profile (figure 2, lower panel) shows significant inner shell absorption 
components to the low energy wing of the ls-2p resonance line which explains why 
simply identifying the absorption near 7 keV with the 6.7 keV rest energy of Fe 
XXV gives too low a velocity. 

Although individually weaker than the Fe absorption, the combination of res¬ 
onance lines of He- and H-like Mg, Si, S and Ar in the broadband spectral fit 
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observed energy (keV) 

Fig. 1 Ratio of EPIC pn data to a simple power law continuum for the 2001 XMM-Newton 
observation of PG121H-143 showing a deep absorption line near 7 keV and additional structure 
between and 4 keV. Deriving an outflow velocity requires the correct identification of 
the individual absorption lines, which ideally requires spectral modelling with a photoionised 
absorber 


is evidently driving the spectral fit. That conclusion is confirmed with Gaussian 
hts to corresponding absorption features in figure 1, which hnd a weighted ob¬ 
served blueshift of 0.055 ± 0.05 and outflow velocity (at the AGN redshift) of 
t; ~ 0.14 ± 0.01c, a value consistent with that found from spectral modelling, but 
significantly higher than from simply identifying the ~7 keV absorption line with 
the resonance ls-2p transitions of either FeXXV or FeXXVI. 

An interesting by-product of the XSTAR modelling in the above case is that 
the observed broadening of the ~ 7 keV absorption line does not require high tur¬ 
bulence (we used grid 25 with vturb of 200 km s“^) or an accelerating/decelerating 
flow. Instead, intrinsically narrow absorption components remain consistent with 
a radial outflow, coasting post-launch. 


2.2 Mass rate and mechanical energy in the PG1211+143 outflow 

Although the detection of high-speed winds in a substantial fraction of bright AGN 
suggests most such flows have a large covering factor, PG1211-I-143 is one of very 
few where a wide angle flow has been demonstrated directly. 

Using stacked data from 4 XMM-Newton observations between 2001 and 2007, 
Pounds & Reeves (2007, 2009) examined the relative strength of ionized emission 
and absorption spectra modelled by XSTAR to estimate the covering factor and 
collimation of the outflowing ionized gas. The summed pn data of PG1211+143 
also shows a well defined P Cygni prohle in the Fe K band (TBTF 3), the classical 
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Fig. 2 (top) A Gaussian fit to the ~7 keV absorption feature finds a line energy of 7.06di0.02 
keV with (Ic) width lOOibSO eV. Identification with the Fe XXV ls-2p resonance line (6.70 keV 
rest energy) gives an outflow velocity v 0.12 ± 0.01c. (lower) Alternative modelling with a 
photoionised gas over the wider 1-10 keV spectral band yields a good fit with a relatively high 
column density A^^'^3.2di0.7xl0^^ cm~^, moderate ionisation parameter log^ = 2.7di0.1 erg 
cm s ~^5 and outflow velocity of v~0.15di0.01c. The Fe XXV absorption line profile is seen to 
include lower energy components due to the addition of one or more L-shell electrons, showing 
why the simple Gaussian fit gives too low a velocity 
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Fig. 3 The PCygni profile of Fe XXV from stacked XMM-Newton pn observations of 
PG1211+143 is characteristic of a wide angle outflow. The comparable equivalent width of 
emission and blue-shifted absorption components indicates the highly ionized outflow has a 
large covering factor. From Pounds and Reeves 2009 


signature of an outflow, with emission and absorption components of compara¬ 
ble equivalent width. Both methods indicated a covering factor 6(= 17/27r) of 
0.75±0.25. Analysis of a Suzaku observation of PG1211-1-143 gives a similar result 
(Reeves et al. 200TBTF), with an intrinsic emission component of ~ 6.5 keV and 
width of (T ~ 250 eV, corresponding to a flow cone of half angle ~ 50 deg, assuming 
velocity broadening in a radial flow. 

The outflow mass rate and mechanical energy can then be estimated, since for 
a uniform radial outflow of velocity v the mass rate is: 


Mout — ATvbnr^mpV, (9) 

where n is the gas density at a radial distance r, and mp is the proton mass. 

The observed values for PG1211-I-143 And a mass loss rate of Mout ~ 7 x 10^^ 
gm s“^ (~ 2.5M0 yr“^), and mechanical energy ~ 4.5 x 10^^ erg s“^ (Pounds & 
Reeves 2009). 

The mass loss rate is comparable to the Eddington accretion rate IW^Edd = 
1.3M0 yr“^ for a supermassive black hole of mass ~ 4 x 1O^M0 accreting at 
an efficiency of 10%, while the outflow mechanical energy is only ~ 6% of the 
Eddington luminosity, close to that predicted by continuum driving (equation 5 
in Section 3 below). As noted elsewhere that energy flow rate would be more than 
sufficient to unbind the gas of the host galaxy bulge if all its energy were efficiently 
communicated. 
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2.3 High speed winds are common 


The evidence for high velocity winds as an important property of AGN remained 
dependent on the prototype case of PG1211+143 for several years, with fast out¬ 
flows in two BAL AGN (Chartas et al. 2002) and in the most luminous low redshift 
QSO PDS 456 (Reeves et al. , 2003 O’Brien et al. 2005) seen as rare objects. That 
began to change with the detection of a highly significant outflow of velocity ~ 0.1c 
in the Seyfert 1 galaxy IC4329A (Markowitz et al. 2006), and several outflow de¬ 
tections in the range ~ 0.14 — 0.2c in multiple observations of Mrk 509 (Dadina 
et al. 2005). A review in 2006 (Cappi et al. 2006) listed 7 non-BAL objects with 
outflows of v~ 0.1c and several with red-shifted absorption lines. 

A major step forward came with the results of an XMM-Newton archival search 
of bright AGN by Tombesi et al. (2010), hnding strong statistical evidence in 15 of 
42 radio-quiet objects of blue-shifted iron K absorption lines, identification with 
FeXXV or XXVI resonance absorption lines implying ultrafast outflow (UFO) 
velocities up to ~0.3c, and clustering near v ~ 0.1c. A later analysis based on 
broad-band modelling with XSTAR photoionized grids (Tombesi et al. 2011) led 
to several revised velocities and conhrmed that the outflows were typically highly 
ionized, with log^ ~ 3 — 6 erg cms“^, with column densities in the range ~ 
10^2 _ ]^q 24 similar search of the Suzaku data archive (Gofford et al. 2013) 

yielded a further group of UFO detections, hnding signihcant absorption in the Fe 
K band in 20 (of 51) AGN with velocities up to ~ 0.3c and a hatter distribution 
than the XMM-Newton sample, but a median value again v~ 0.1c. 

Figure 4 brings together the results from the spectral modelling analyses of 
the XMM-Newton and Suzaku surveys. We follow Tombesi et al. (2011) in dehning 
UFOs as having outhow velocities greater than 10"^ km s“^, to discriminate against 
WAs or post-shock hows (Sections 2.4 and 4). The velocity plot shows a peak 
at ~ 0.1c, with a tail extending to ~ 0.3c. In terms of the continuum-driving 
Black Hole Winds model (King & Pounds 2003) the higher velocities would imply 
a higher value of the accretion efhciency rj, with the future potential for such 
observations to provide a measure of black hole spin. Equation (22) also suggests 
the low velocity tail in both the Tombesi et al. and Gofford et al. distributions 
could relate to primary outhows formed at a higher accretion ratio (but see Section 
3.1). 

Figure 4 also shows the distribution of ionization parameter and absorption 
column density from the surveys of Tombesi et al. (2011) and Gofford et al. 
(2013). The high ionization parameter, peaking near log^ ~ 4, explains why the 
detection of UFOs has been almost exclusively from X-ray observations in the Fe 
K band, leaving open the possibility that fully ionized outflows (also consistent 
with continuum driving) will become detectable when the AGN luminosity (and 
hence ionization) falls. In assessing observational data it is important to note that 
for a radial outflow the observed column density is a line-of-sight integration 
over the flow time, dominated by the higher density at small radii, while the 
ionization parameter is governed by the current AGN luminosity. The column 
density, which generally lies below Njj ~ 10^"^ cm“^, can vary rapidly and turns 
out to be a powerful diagnostic of the flow history and dynamics. We return to 
the observability of UFOs in Section 3.3. 
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Fig. 4 Distribution of outflow velocities, ionization parameter (erg cm s and column den¬ 
sity (cm“^) obtained from modelling the individual spectra from extended observations of 
type 1 AGN in the XMM-Newton and Suzaku data archives (Tombesi et al. 2011, Gofford et 
al. 2013). The red lined histogram refers to lower limits in column density 


2.4 Evidence for a shocked flow 

The mechanical energy in a fast wind, such as that in PG1211+143, was noted in 
Section 2.2 to be incompatible with the continued growth of the black hole and 
stellar bulge of the host galaxy, unless the flow is short-lived or the coupling of 
wind energy to bulge gas is highly inefhcient. A recent XMM-Newton observation 
of the narrow-line Seyfert galaxy NGC 4051 has provided the first evidence of a 
fast ionized wind being shocked, with subsequent strong cooling leading to most 
of the initial flow energy being lost before it can be communicated to the bulge 
gas. We outline a possible scenario for that event below. 

NGG 4051 was found in the XMM-Newton archival search to have a high velocity 
wind during an observation in 2002 when the source was in an unusually low state, 
the initial identification with Fe XXVI Lyman-a in Tombesi et al. (2010) indicating 
a velocity of ~0.15c. In a full spectral fit (Tombesi et al. 2011) identification with 
Fe XXV was preferred, with an increased velocity ~0.20c. Significantly, in a 2001 
observation of NGC 4051, when the X-ray flux was much higher, a strong outflow 
was detected at ~6000 km s“^, but no ultra-fast wind was seen. 
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Fig. 5 The outflow velocity and ionization parameters for 6 XSTAR photoionised absorbers 
used to fit the RGS and EPIC spectra of NGC 4051, together with a high point representative 
of the pre-shock wind, show the linear correlation expected for a mass-conserved cooling flow 
(see Pounds and King 2013) 


It seems that the detection of a UFO in NGC 4051 is unusually dependent on 
the source flux, with evidence for a high velocity wind {v ~ 0.13c) again found 
only during periods when the ionizing continuum was low during a further XMM- 
Newton observation in 2009 (Pounds & Vaughan 2012). An additional factor may 
be the low redshift (z=0.00234) of NGC 4051, which makes a high velocity wind 
more difficult to detect with current observing facilities. 

The 600 ks XMM-Newton observation of the Seyfert 1 galaxy NGC 4051 in 
2009, extending over 6 weeks and 15 spacecraft orbits, broke new ground by find¬ 
ing an unusually rich absorption spectrum with multiple outflow velocities, in both 
RGS (den Herder et al. 2001) and EPIC spectra, up to ~ 9000 kms~^ (Pounds 
& Vaughan 2011a). Inter-orbit variability is seen in both absorption and emission 
lines, with strong recombination continua (RRC) and velocity-broadened reso¬ 
nance lines providing constraints on the dynamics and geometry of the putative 
post-shock flow (Pounds and Vaughan 2011b, 2012). 


2.5 A self-consistent model for the shocked wind in NGC 4051 

More complete modelling of both RGS and EPIC pn absorption spectra of NGC 
4051 found a highly significant correlation of outflow velocity and ionization state 
(figure 5), as expected from mass conservation in a post-shock flow (King 2010, 
Pounds & King 2013). The additional analysis also found a range of column den¬ 
sities to be required by the individual XSTAR absorption components, suggesting 
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Fig. 6 Fe K profiles from observations of NGC 4051 several days apart show an increased 
level of ionization coinciding with a hard X-ray flare (data from Pounds and Vaughan 2012). 
The ratio of resonance absorption lines of Fe XXV and Fe XXVI is a sensitive measure of the 
ionization state of the absorbing gas 


an inhomogeneous shocked flow, perhaps with lower ionization gas clumps or fila¬ 
ments embedded in a more extended, lower density and more highly ionized flow. 

Theoretical considerations suggested a key factor in determining the structure 
of the post-shock flow was likely to be the cooling time, as discussed in more 
detail in Section 4. In particular, the fate of a fast wind depends on the distance 
it travels before colliding with the ISM or slower-moving ejecta, with Compton 
cooling dominating for a shock occurring sufficiently close to the AGN continuum 
source. 

Importantly, flux-linked variations in the ratio of FeXXV to Fe XXVI absorp¬ 
tion in the 2009 XMM-Newton observation (figure 6) provided a measure of the 
Compton cooling time, the mean flow speed then determining the shell thick¬ 
ness of the hotter, more highly ionized flow component. The detection of strong 
recombination continua (RRC) in the soft X-ray spectra furthermore suggested 
an increasing density in the decelerating post-shock flow, with two-body cooling 
becoming increasingly important. 

To pursue that idea we note that at the (adiabatic) shock the free-free (thermal 
bremsstrahlung) and Compton cooling times are 

tff~3x s = 20-^ yr (10) 

N Mrm ^ ’ 


and 


to = 10"^ 


Ms 


yr 


( 11 ) 


respectively (see King et al. 2011: here T,N are the postshock temperature and 
number density, Rie is the shock radius in units of 10^® cm, M^ is the black hole 
mass in units of IO^Mq, and m ~ 1 is the Eddington factor of the mass outflow 
rate). 

After the adiabatic shock, the gas cools rapidly via inverse Compton cooling, 
while its density rises as A oc T~^ (pressure is almost constant in an isothermal 
shock), and 


tff oc 


N 


( 12 ) 
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which means that the free-free cooling time decreases sharply while the Compton 
time does not change. Eventually free-free (and other atomic two-body processes) 
become faster than Compton when T has decreased sufficiently below the original 
shock temperature Ts ~ 1.6 x 10^° K. From (11011111) above this requires 



(13) 


or 


T < 2 X 10^ K. 


(14) 


The temperature of ionization species forming around a few keV is therefore 
likely to be determined by atomic cooling processes rather than Compton cool¬ 
ing. The strong recombination continua in NGC 4051 (Pounds & Vaughan 2011b, 
Pounds & King 2013) are direct evidence for that additional cooling, with the RRC 
flux yielding an emission measure for the related flow component. In particular, 
the onset of strong two-body cooling results in the lower-ionization, lower-velocity 
gas being conhned in a relatively narrow region in the later stages of the post¬ 
shock flow. The structure and scale of both high and low ionization flow regions 
can be deduced from the observations and modelling parameters. 

For the highly ionized post-shock flow, the iron Ly-a to He-a ratio will be 
governed by the ionizing continuum and recombination time. Signihcant varia¬ 
tions in this ratio are found on inter-orbit timescales (Pounds & Vaughan 2012), 
with an example shown in hgure 6. For a mean temperature of ~1 keV, and 
recombination coefficient of 4.6 x 10“^^ cm^s“^ (Verner & Ferland 1996), the 
observed recombination timescale of ~ 2 x 10® s corresponds to an average par¬ 
ticle density of ~ 4 x 10®cm“®. Comparison with a relevant absorption column 
Nfi ~ 4 X 10^^ cm“^ from the XSTAR modelling indicates a column length scale of 
~ 10^® cm. Assuming a mean velocity of the highly ionized post-shock flow of 6000 
km s“^, the observed absorption length corresponds to a flow time ~ 1.7 x 10^ s 
(0.6 yr). Equation (11) finds a comparable cooling time for NGC 4051 at a shock 
radius R ~ 10^^ cm. 

For the low-ionization flow component, decay of strong RRC of NVII, CVI 
and CV (Pounds & Vaughan 2011b, Pounds & King 2013), occurs over ~ 2 — 6 
days. With an electron temperature from the mean RRC profile of ~ 5 eV, and 
recombination coefficient for CVI of ~ 10“^^ cm®s“^ (Verner & Ferland 1996), 
the observed RRC decay timescale corresponds to a (minimum) electron density 
of ~2 X 10® cm“®. A column density of 1.5 x 10^^ cm“^ from modelling absorption 
in the main low-ionization flow component then corresponds to an absorbing path 
length of 7 X 10^'^ cm. 

The RRC emission flux provides a consistency check on the above scaling. 
Assuming solar abundances, and 30 percent of recombinations direct to the ground 
state, a CVI RRC flux of ~ 10“® photons cm“^ s“^ corresponds to an emission 
measure of ~ 2 x 10®^ cm“®, assuming a Tully-Fisher distance to NGC 4051 of 
15.2 Mpc. With a mean particle density of ~ 2 x 10® cm“® the emission volume 
[A'kR?AR) is ~ 5 X 10^® cm® . Assuming a spherical shell geometry of the flow, 
with fractional solid angle b, shell thickness AR ~ 7 x 10®'^ cm, and shell radius 
R ~ 10^^ cm, the measured RRC flux is reproduced for b ~ 0.5. 

Although this excellent agreement may be fortuitous given the approximate 
nature and averaging of several observed and modelled parameters, the mutual 
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shock width ~ 10^^ cm 

Fig. 7 Schematic view of the shock pattern resulting from the impact of a black hole wind 
(blue) on the interstellar gas (red) of the host galaxy. The accreting supermassive black hole 
drives a fast wind (velocity v ~ riclm 0.1c), whose ionization state makes it observable in 
X-ray absorption lines. It collides with the ambient gas in the host galaxy and is slowed in a 
strong shock. The inverse Compton effect from the quasar’s radiation field rapidly cools the 
shocked gas, removing its thermal energy and strongly compressing and slowing it over a very 
narrow radial extent. In the most compressed gas, two—body cooling becomes important, and 
the flow rapidly cools and slows over an even narrower region. In NGC 4051 this region is 
detected in the soft X-ray spectrum, where absorption (and emission) are dominated by the 
lighter metals. The cooled gas exerts the preshock ram pressure on the galaxy’s interstellar 
gas and sweeps it up into a dense shell (‘snowplow’). The shell’s motion then drives a milder 
outward shock into the ambient interstellar medium. This shock ultimately stalls unless the 
SMBH mass has reached the value Mo- satisfying the M — a relation (from Pounds and King 
2013) 


consistency of absorption and emission of the photoionized flow is encouraging. 
Given that only blueshifted RRC emission is seen, b ~0.5 is consistent with a 
wide-angle flow, visible only on the near side of the accretion disc. 

Figure 7 illustrates the main features of the overall NGC 4051 outflow, a fast 
primary wind being shocked at a radial distance of order O.lpc, within the zone 
of influence of an SMBH of 1.7 x 10®Mq. The initially hot gas then cools in the 
strong radiation held of the AGN, with a Compton cooling length determining the 
absorption columns of Fe and the other heavy metal ions. Two-body recombination 
provides additional cooling as the density rises downstream, eventually becoming 
dominant. Absorption (and emission) in the soft X-ray band is located primarily 
in this thinner outer layer of the post-shock flow. 
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It is interesting to note that similar shocking of fast outflows provides a natural 
link between UFOs and the equally common ‘warm absorbers’ in AGN (Tombesi et 
al. 2013). While the onset of strong two-body cooling, resulting in the intermediate 
column densities being small, might explain why evidence for intermediate-flow 
velocities has awaited an unusually long observation of a low mass AGN, the 
accumulated ‘debris’ of shocked wind and ISM could be a major component of the 
‘warm absorber’. See Section 7.4.1 for a discussion. 


2.6 Variability of UFOs 

While it is likely that powerful winds blow continuously in AGN in rapid growth 
phases, it is important to note that the existing observations of UFOs are restricted 
to bright, low-redshift AGN, z < 0.1, where the X-ray fluxes are sufficient to yield 
high quality spectra. Repeated observation of several bright AGN frequently show 
changes in the equivalent width in the primary Fe K absorption line. 

Variability in the strength of blueshifted Fe-K absorption over several years in 
PG1211-1-143 was first noted in a comparison of the initial XMM-Newton and Chan¬ 
dra observations (Reeves et al. 2008), and conhrmed by repeated XMM-Newton 
observations (Pounds and Reeves 2009). Multiple observations of the luminous 
Seyfert 1 galaxy Mrk 509 (Gappi et al. 2009) found variations in both inten¬ 
sity and blueshift of Fe K absorption lines. The archival searches provide the 
most comprehensive variability data, with repeated observations of several AGN 
demonstrating that variability of absorption line equivalent width (EW) over sev¬ 
eral years is common. More rapid variability in EW, over a few months, is reported 
in the XMM-Newton archive for Mrk 509, Mrk 79 and Mrk 841, with both velocity 
and EW change in <2 days for Mrk766. 

In addition, the ‘hit rate’ of UFOs for multiply-observed AGN in the archival 
XMM-Newton data search (Tombesi et al. 2010) was relatively low, being 1 of 6 
observations for NGG4151, MCG-6-30-15 (0/5), Mrk509 (3/5), NGC 4051(1/2), 
Mrk79 (1/3), Mrk205 (1/3), and Mrk290 (1/4). Overall, though 101 suitably ex¬ 
tended observations yielded 36 narrow absorption line detections in the Fe K band, 
only 22 were observed at > 7 keV. While the UFO ‘hit rate’ of ~ 22% is a lower 
limit set by the sensitivity of available exposures it seems clear that the fast out¬ 
flows currently being detected in low-redshift AGN are far from continuous. 


3 BLACK HOLE WINDS 

3.1 The Eddington Accretion Ratio in AGN 

We have seen from Section 2 that a large fraction of observed AGN show ultrafast 
outflows. Unless we view every AGN from a very particular angle (so implying a 
much larger total population) this must mean that these winds have large solid an¬ 
gles 47r6 with fe ~ 1, i.e. they are quasispherical. We recall that UFOs are observed 
to have total scalar momenta ~ L^dd/c, where LEdd is the Eddington luminosity 
of the SMBH. We can argue that on quite general grounds, SMBH mass growth is 
likely to occur at accretion rates close to the value MEdd = L^dd/v^^ which would 
produce this luminosity. As we noted earlier, the Soltan (1982) relation shows 
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that the largest SMBH gained most of their mass by luminous accretion, i.e. dur¬ 
ing AGN phases. But the fraction of AGN among all galaxies is small, strongly 
suggesting that when SMBH grow, they are likely to do so as fast as possible. The 
maximum possible rate of accretion from a galaxy bulge with velocity dispersion 
a is the dynamical value 

Mdyn ^ (15) 

where fg is the gas fraction. This rate applies when gas which was previously in 
gravitational equilibrium is disturbed and falls freely, since one can estimate that 
the gas mass was roughly Mg ~ a'^fgR/G. Once this is destabilized it must fall 
inwards on a dynamical timescale tdyn ~ R/u. This gives the result m , since 

Mdyn ~ Mg/tdyn- 

Numerically we have 

Mdyn — 280 cr|oo M0yr ^ ( 16 ) 


where we have taken fg 
have 


0.16, the cosmological baryon fraction of all matter. We 


MEdd = 


^Edd 

rjc^ 


AttGM 

uric 


(17) 


where I/Edd is the Eddington luminosity and k is the electron scattering opacity. 
With t ; = 0.1 and black hole masses M close to the observed M — a relation © 
we find 

MEdd — 4.4 cr2ooM0yr ^ (18) 

and an Eddington accretion ratio 


Mdyn _ _ 54 

MEdd ~ <^200 “ ’ 


(19) 


Thus even dynamical infall cannot produce extremely super—Eddington accre¬ 
tion rates on to supermassive black holes. But the rate m is already a generous 
overestimate, since it assumes that the infalling gas instantly loses all its angular 
momentum. Keeping even a tiny fraction of this instead forces the gas to orbit the 
black hole and form an accretion disc, which slows things down drastically. Gas 
spirals inwards through a disc on the viscous timescale 


1 

a 



( 20 ) 


where a ~ 0.1 is the Shakura-Sunyaev viscosity parameter, while the disc aspect 
ratio H/R is almost constant with radius, and typically close to 10“^ for an AGN 
accretion disc (e.g. Collin-Souffrin & Dumont 1990). Then Gisc approaches a Hub¬ 
ble time even for disc radii of only 1 pc. Gas must evidently be rather accurately 
channelled towards the SMBH in order to accrete at all, constituting a major 
problem for theories of how AGN are fuelled. 

Given all this, we see that while it is possible for AGN accretion rates to reach 
Eddington ratios m ~ 1, significantly larger ones are unlikely unless the SMBH 
mass is far below the M — a value appropriate to the galaxy bulge it inhabits. In 
other words, only relatively modest values m ~ 1 of the Eddington ratio are likely 
in SMBH growth episodes. 
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Indirect evidence supporting this view comes from stellar-mass compact binary 
systems. The dynamical rate is relatively much larger here, as the equivalent of 
(HB is ^ ^ Ivl‘ 2 !P with Uorb the orbital velocity of a companion star in 

a binary of period P. This implies rates approaching a solar mass per few hours in 
many cases, if dynamical accretion ever occurs. These systems have highly super- 
Eddington apparent luminosities, probably as the result of geometric collimation 
(cf King et al. 2001). But significantly there are no obvious AGN analogues of the 
ultra-luminous X-ray sources (ULXs), suggesting that Eddington ratios m ^ 1 are 
very unusual or absent in AGN. 


3.2 Eddington Winds 

Given this, we can crudely model the UFOs discussed in Section 2 as quasi- 
spherical winds from SMBH accreting at modest Eddington ratios m = M/M-^dd ~ 
1. Winds like this have electron scattering optical depth r ~ 1, measured inwards 
from infinity to a distance of order the Schwarzschild radius (cf eq. [27| below). So 
on average every photon emitted by the AGN scatters about once before escap¬ 
ing to infinity. Since electron scattering is front—back symmetric, each photon on 
average gives up all its momentum to the wind, and so the total (scalar) wind 
momentum should be of order the photon momentum, or 

( 21 ) 

c 

where v is the wind’s terminal velocity. The winds of hot stars obey relations like 
this. For accretion from a disc, as here, the classic paper of Shakura & Sunyaev 
(1973) finds a similar result at super-Eddington mass inflow rates: the excess 
accretion is expelled from the disc in a quasispherical wind. 

Equation (nm now directly gives the wind terminal velocity as 

n~^c~0.1c. (22) 

m 

From eq. (|22p we get the instantaneous wind mechanical luminosity as 

isHwind = ~ ~ ^LEdd O.OSLEdd. (23) 

2 c 2 2 

This relation turns out to be highly significant (see Sections 5.3, 7.3). 

Ohsuga & Mineshige (2011) show in detail that winds with these properties 
(their Models A and B) are a natural outcome of mildly super-Eddington accre¬ 
tion. In particular their Model A and B winds are predicted (cf their Figure 3) to 
have mechanical luminosities ~ O.lZ/Edd, in rough agreement with equation (1231) . 

Compared with the original disparity Ebh = rjMc^ ~ 2000i5gas between black 
hole and bulge gas binding energies outlined in the Introduction, we now have a 
relation ^ 

EsHwind — ~ lOOEgas (24) 

between the available black hole wind mechanical energy and the bulge binding en¬ 
ergy. Although the mismatch is less severe, it still strongly suggests that the bulge 
gas would be massively disrupted if it experienced the full mechanical luminosity 
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emitted by the black hole for a significant time. So the coupling of mechanical en¬ 
ergy to the host ISM cannot be efficient all the time (see the discussion in Section 
2). We show how this works in Section 4. 


3.3 Observability 


As we have seen, observations frequently give the hydrogen column density Nu 
through a UFO wind from the X-ray absorption spectrum. We can show that this 
quantity determines whether a given UFO wind is observable or not. Using m 
in the mass conservation equation 

Mw = 47rfer^wp(r), (25) 


where p(r) is the mass density, we find the equivalent hydrogen column density of 
the wind as 


Nh 



Mw 

Airr'^bv 


dr = 


-bEdd 

AnbrnpRi^cv^ 


(26) 


where i?in is its inner radius, mp is the proton mass, and we have used (ED at the 
last step. From the definition of L^dd we find the wind electron scattering optical 
depth 


T = NpjfJJ' 


GM 

bv'^Rin 


(27) 


with ax — ump the Thomson cross-section. This shows self-consistently that the 
scattering optical depth r of a continuous wind is ~ 1 (cf King & Pounds 2003, 
equation 4) at the launch radius i?iaunch — GM/bv^ = {<?/ 2 bv^)Rs ~ 50i?s. 

The measured values of Npj (Tombesi et al. 2011, Gofford et al. 2013, Figure 
4) are always smaller than the value Npj ~ \/<jt — 10^"^ cm“^ for a continuous 
wind, and actually lie in the range N 22 ~ 0.3 — 30, where N 22 = cm“^. 

It is perhaps not surprising that observations do not show any UFO systems with 
Npi > 10^“^ cm“^. These AGN would be obscured at all photon energies by electron 
scattering, and perhaps difficult to see. Although such systems might be common, 
we probably cannot detect them. To have a good chance of seeing a UFO system 
we need a smaller N}j, so from ED the inner surface 7?in of the wind must be 
larger than Aiaunch- This is only possible if all observed UFOs are episodic, i.e. we 
see them some time after the wind from the SMBH has switched off. In this sense 
UFOs are more like a series of sporadically-launched quasispherical shells than a 
continuous outflow. The Nu value of each shell is dominated by the gas near its 
inner edge (cf (I26II , so we probably at most detect only the inner edge of the most 
recently-launched shell. We can quantify this by setting i?in = wtoff, where toff is 
the time since the launching of the most recent wind episode ended. Using ED 
gives 


GM 


3Mr 

bv^-^^N22 


months. 


(28) 


where wo.i = n/O.lc. 

As seen in Figure 4, all UFOs have N22 ~ 0.3 — 30, and most of the SMBH 
masses are ~ IO^Mq. Evidently the launches of most observed UFO winds halted 
weeks or months before the observation. At first sight this is surprising. The 
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strength of the characteristic blueshifted absorption features defining UFOs is 
closely related to Njj. These features would be still stronger if there were UFOs 
with N 22 > 100, but none are seen. We note from (I28II that observing a UFO like 
this would require us to catch it within days of launch. Given the relatively sparse 
coverage of X-ray observations of AGN this is unlikely. So the apparent upper 
limit to the observed may simply reflect a lack of observational coverage, and 
implies that most UFOs are short-lived 

The lower limit to in the Tombesi et al. sample is also interesting. Once 
Nfj is smaller than some critical value, any blueshifted absorption lines must 
become too weak to detect. The strongest are the resonance lines of hydrogen- 
and helium-like iron, which have absorption cross-sections crpe ~ 10“^® cm^. Given 
the abundance by number of iron as Zpe = 4 x 10“® times that of hydrogen, the 
condition that one of these lines should have significant optical depth translates 
to ZpeNno'Fe > 1 or N 22 > 2.5. This is similar to the lowest observed values. From 
(I28II this means that current observations cannot detect UFO winds launched 
more than a few months in the past, because the blueshifted iron lines will be 
too weak. Even these observed UFOs should gradually decrease their Njj and 
become unobservable if followed for a few years. We see in Section 5 that the 
UFO wind typically travels ~ IOM 7 pc or more before colliding with the host 
galaxy’s interstellar gas, which takes tcoii/v ~ yr. Finally, a UFO may 

be unobservable simply because it is too strongly ionized, so that no significant 
Nfj can be detected. 

All this means that the state of the AGN seen in a UFO detection does not 
necessarily give a good idea of the conditions required to launch it. In particular, 
the AGN may be observed at a sub-Eddington luminosity, even though one might 
expect luminosities ~ TEdd fo be needed for launching the UFO. This may be the 
reason why AGN showing other signs of super-Eddington phenomena (e.g. narrow¬ 
line Seyfert 2 galaxies) are nevertheless seen to have sub-Eddington luminosities 
most of the time (e.g. NGC 4051; Denney et al. 2009): the rather short wind 
episodes are launched in very brief phases in which accretion is slightly super- 
Eddington, whereas the long-term average rate of mass gain may be significantly 
sub-Eddington. 

In summary, it is likely that current UFO coverage is remarkably sparse. We 
cannot see a continuous wind at all. We can only see an episodic wind shell shortly 
after launch, and then only for a tiny fraction toff/fcoii ~ of ifs 

~ 300 — 3000 yr journey to collision with the host ISM. So it seems that the vast 
majority of UFO wind episodes remain undetected: more AGN must produce them 
than we observe, and the known UFO sources may have far more episodes than 
we detect. 

All this has important consequences for how we interpret observations in dis¬ 
cussing feedback. The most serious is that the most powerful form of feedback - 
from AGN at the Eddington limit producing continuous winds - is probably not 
directly observable at all. 


3.4 The Wind Ionization State, and BAL QSOs 




NR"^ 


The ionization parameter 


(29) 
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essentially fixes the ionization state of a black hole wind wind, and so determines 
which spectral lines are observed. Here Li = liL-^dd is the ionizing luminosity, with 
< 1 a dimensionless parameter specified by the quasar spectrum, and N = pj[imp 
is the number density of the UFO gas. We use dll [251) to get 

^ = 3 X 10‘‘r7o.i^2m“^ = 3 x 10‘^Uo.i^2 (30) 

where = 1^/10“^, and 770.1 = 77 / 0 . 1 . 

This relation shows how the wind momentum and mass rates determine its 
ionization parameter and so its line spectrum as well as its speed v. Given a quasar 
spectrum Lv , the ionization state has to arrange that the threshold photon energy 
defining and the corresponding ionization parameter together satisfy (|30ll. 
This shows that the excitation must be high: a low threshold photon energy (say in 
the infrared) would imply a large value of I 2 , but then (1301) gives a high value of ^ 
and so predicts the presence of very highly ionized species, physically incompatible 
with such low excitation. 

For suitably chosen continuum spectra (PPj) has a range of solutions. A given 
spectrum might in principle allow more than one solution, the applicable one being 
specified by initial conditions. For a typical quasar spectrum, an obvious self- 
consistent solution of (|30|) is I 2 — 1, 711 ~ 1, ^ ~ 3 X 10^. Here the quasar radiates 
the Eddington luminosity. We can also consider situations where the quasar’s 
luminosity has decreased after an Eddington episode but the wind is still flowing, 
with 771 ~ 1. Then the ionizing luminosity 10“^l2TEdd in (1301) is smaller, implying 
a lower /. As an example, an AGN of luminosity 0.3Z/Edd would have ^ ~ 10^. 
This gives a photon energy threshold appropriate to FeXXV and Fe XXVI (i.e. 
^T^threshoid ~ 9 keV). We conclude that Eddington winds from AGN are likely to 
have velocities ~ 0.1c, and show the presence of helium- or hydrogen-like iron in 
accord with the absorption reported in Section 2. Zubovas & King (2013) show 
that this probably holds even for AGN which are significantly sub-Eddington. 

We can see from (1221) that a larger Eddington factor rh is likely to produce 
a slower wind. From comparison with ULXs (see Section 3.1) we also expect the 
AGN radiation to be beamed away from a large fraction of the UFO, which should 
therefore be be less ionized, and as a result more easily detectable than the small 
fraction receiving the beamed radiation. These properties - slower, less ionized 
winds - characterize BAL QSO outflows, perhaps suggesting that systems with 
larger rh > 1 appear as BAL QSOs. Zubovas & King (2013) tentatively confirm 
this idea. 


4 THE WIND SHOCK 

4.1 Momentum- and Energy-Driven Flows 

So far we have only studied the black hole wind. But we know that this wind 
must have a significant effect on the host galaxy when it impacts directly on its 
interstellar medium (ISM). In this Section we model the wind and host ISM as 
roughly spherically symmetric, and consider the effects of deviations from this 
simple picture later. 

The pattern of the wind-ISM interaction (Figure 7) is qualitatively identical to 
that of a stellar wind hitting the interstellar medium around it (see e.g. Dyson & 
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/ rapidly lost to cooling, 

contact discontinuity only ram pressure 

oMBH ('momentum and conserved) 



ISM shock 
(energy conserved) 


Fig. 8 Schematic picture of momentum-driven (top) and energy—driven (bottom) outflows. 
In both cases a fast wind (velocity ~ 0.1c) impacts the interstellar gas of the host galaxy, 
producing an inner reverse shock slowing the wind, and an outer forward shock accelerating 
the swept-up gas. In the momentum—driven case (top), corresponding to the UFOs discussed 
in Section 2, the shocks are very narrow and rapidly cool to become effectively isothermal. 
Only the ram pressure is communicated to the outflow, leading to very low kinetic energy 
(cr/c)I/Edd- Iri an energy-driven outflow (bottom), the shocked regions are much wider and 
do not cool. They expand adiabatically, communicating most of the kinetic energy of the wind 
to the outflow (in simple cases divided in a ratio of about 1:2 between the shocked wind and 
the swept-up gas). The outflow radial momentum flux is therefore greater than that of the 
wind. Momentum-driven conditions hold for shocks confined to within ~ 1 kpc of the AGN, 
and establish the M — a relation gg (King, 2003; King, 2005). Once the supermassive black 
hole mass attains the critical M — a value, the shocks move further from the AGN and the 
outflow becomes energy—driven. This produces the observed large-scale molecular outflows 
which probably sweep the galaxy clear of gas. (From Zubovas & King, 2012a). 

Williams 1997). The black hole wind (shown in bine) is abruptly slowed in an inner 
(reverse) shock where the temperature approaches ~ 10^^ K if ions and electrons 
reach equipartition (but see the discussion below). The shocked wind gas acts like a 
piston, sweeping up the host ISM (shown in red) at a contact discontinuity moving 
ahead of it. Because this swept-up gas moves supersonically into the ambient ISM, 
it drives an outer (forward) shock into it (see Figs. 7 and 8 [top]). 

The dominant interaction here is the reverse shock slowing the black hole wind, 
and injecting energy into the host ISM. The nature of this shock differs sharply 
depending on whether or not some form of cooling (typically radiation) removes 
signihcant energy from the hot shocked gas on a timescale shorter than its flow 
time. If the cooling is strong in this sense (‘momentum-driven flow’), most of the 
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preshock kinetic energy is lost (usually to radiation). The very rapid cooling means 
that the shocked wind gas is highly compressed, making the postshock region 
geometrically narrow (see the upper part of Figure 8). This kind of narrow, strongly 
cooling region ia often idealised as a discontinuity, know as an ‘isothermal shock’ 
(cf Dyson & Williams, 1997). As momentum must be conserved, the postshock gas 
transmits just its ram pressure (ED to the host ISM. We will see that this amounts 
to transfer of only a fraction ~ cr/c ~ 10“^ of the mechanical luminosity i?BH wind — 
O.OSZ/Edd (cf equation I23II to the ISM. In other words, in the momentum-driven 
limit, only energy 

2 

Smom ~ ^F^BHwind ~ ^ ~ ^ X ~ O.lEgas (31) 

is injected into the bulge ISM, i.e. about 10% of the bulge gas binding energy 
fgMi,a^ for black holes close to the M — a relation (there is now an implicit factor 
0 - 200 /Ms cci the rhs). Thus momentum-driven flows do not threaten the bulge’s 
integrity. Indeed we will see that they never interact with most of it, so there is 
no danger that the black hole will drive away the gas and suppress accretion. A 
momentum-driven regime is a stable environment for black-hole mass growth. 

In the opposite limit where cooling is negligible, the postshock gas retains all 
the mechanical luminosity 

7?wind — O.OSZ/Edd — lOOFgas (32) 

(cf equation 1241) thermalized in the shock, and instead expands adiabatically into 
the ISM. The postshock gas is now geometrically extended (see the lower part of 
Figure 8), unlike the momentum-driven (‘isothermal’) case. This ‘energy-driven 
flow’ is much more violent than momentum-driven flow. The estimate (13211 is for a 
black-hole mass near the M—a relation; a hole with mass a factor of 100 below this 
would already unbind the bulge in doubling its mass. Unless the shock interaction 
is markedly aspherical, a black hole in an energy—driven environment is unlikely 
to reach observed SMBH masses. 

Given these starkly different outcomes we must decide under what conditions 
we have momentum- or energy-driven outflows. Simple estimates immediately 
show that ordinary atomic two-body processes have no signihcant effect in cool¬ 
ing the wind shock. But the wind shock is exposed to the radiation held of an 
Eddington-accreting supermassive black hole. This has a characteristic tempera¬ 
ture of no more than ~ 10^ K, far lower than the wind’s immediate post-shock 
temperature of ~ 10^° — 10^^ K. Electrons in the postshock gas lose energy to 
these photons through the inverse Compton effect (cf Ciotti & Ostriker, 1997), at 
a rate dependent on the radiation density. Eor wind shocks close to the SMBH, 
the accretion radiation held is intense enough that this effect cools the postshock 
wind gas in less than the momentum-driven how time ~ Rja (see below), and we 
are self-consistently in the momentum-driven regime, provided that the postshock 
gas is in equipartition, i.e. that electron and ion temperatures remain effectively 
equal. (We consider this further in Section 4.2 below.) 

For shocks at larger radii R the radiation energy density decreases as 
increasing the cooling time as R?. The how time increases only as A, so for R 
greater than the critical cooling radius 

Rc ~ 500Mg^^CT2OO pc 


(33) 
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(King, 2003, 2005; King et al 2011, Znbovas & King 2012b) the cooling time 
is longer than the flow time, and the flow must be energy-driven. So we have the 
general result that momentum-driven flows are confined to a small region R < Rq-, 
while energy-driven flows must be large-scale. This is just as one would expect, 
given that a momentum—driven flow allows stable black hole mass growth, while 
an energy-driven one is likely to expel most of the bulge gas. 

It is plausible then that the observed UFO winds can lead to momentum¬ 
driving through strongly cooled shocks close to the SMBH. In this picture all but a 
small fraction of the mechanical luminosity (1231) of the black hole wind is eventually 
radiated away as an inverse Compton continuum with characteristic photon energy 
~ 1 keV. Pounds & Vaughan (2011) report a possible detection of this spectral 
component in the Seyfert 1 galaxy NGC 4051. As required for consistency, the 
luminosity of this component is comparable to the expected mechanical luminosity 
(1231) of the wind in that system. 

Cooling shocks are called ‘isothermal’ because the gas temperature rapidly re¬ 
turns to something like its preshock value. Momentum conservation requires that 
the gas is also strongly slowed and compressed as it cools. So the postshock veloc¬ 
ity of the X-ray emitting gas should correlate with its temperature (or roughly, 
ionization) while Compton cooling is dominant. Once this has compressed the gas 
sufficiently, two-body processes such as free-free and bound-free emission must 
begin to dominate, since they go as the square of the density, and their cooling 
times decrease with temperature also (Pounds & King 2013). Section 2.5 above 
shows that there is direct observational evidence for both of these effects in NGC 
4051. So this object (uniquely) shows three signatures of a cooling shock: an in¬ 
verse Compton continuum, an ionization-velocity correlation, and the appearance 
of two-body processes in the spectrum. 


4.2 Shock Cooling 

Cooling (or the lack of it) has a defining effect on the physics of the interaction 
between the black hole wind and the host ISM, so we must check the simple 
picture above. In particular the inverse Compton effect acts only on electrons, but 
the energy of the postshock gas is initially almost all in its ions. We assumed above 
that the electron and ion temperatures quickly come into equipartition after the 
shock, allowing the inverse Compton effect to drain energy from the ions. 

This assumption can be questioned. Faucher-Giguere & Quataert (2012) show 
that if the only process coupling electrons and ions is Coulomb collisions, there is a 
significant parameter space where equipartition does not occur, although they do 
not rule out substantial momentum-conserving phases. An important considera¬ 
tion here is that many processes other than direct Coulomb collisions may rapidly 
equilibrate electron and ion temperatures. Faucher-Giguere & Quataert (2012) 
attempt to put limits on the incidence of such collisionless coupling by appeal¬ 
ing to observations of the solar wind, but this is an area of considerable physical 
uncertainty. 

Another way of using observations to decide if shock cooling is effective is to 
look for the inverse Compton spectral component directly revealing the cooling. 
Bourne et al. 2013) argue that the apparent lack of such a component in most AGN 
spectra rules out cooling shocks. But we recall from Section 3.3 that the coverage 
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of UFOs is extremely sparse. Actually observing a collision and so the inverse 
Compton emission is inevitably a very rare event. It appears that observationally 
ruling out Compton shock cooling is so far inconclusive. 


5 THE M - a RELATION 

5.1 Reaching M — a\ the Momentum-Driven Phase 

We are now equipped to discuss the impact of a UFO on the host interstellar 
gas. We already noted that the wind impact implies a pair of shocks each side of 
the contact discontinuity between the wind and the host ISM. Initially the wind 
shock is close the hole, and we assume that inverse Compton cooling from the 
AGN radiation field cools it rapidly and puts the flow in the momentum-driven 
regime. The region of gas between the wind shock and the contact discontinuity, 
where it impacts and sweeps up the host ISM, is very narrow (cf the upper panel of 
Figure 8). The outer shock accelerating the ISM is also strongly cooled, so that the 
‘snowplow’ region of swept-up ISM is narrow as well. So we can treat the whole 
region between the inner and outer shocks as a single narrow, outward-moving 
gas shell, whose mass grows as it sweeps up the host ISM (see Fig. 9). 

As a simple model of a bulge, we assume that the matter of the host galaxy is 
distributed with an isothermal profile of velocity dispersion a, with mass density 


p(r) = 


so that the mass within radius R is 


2'KGr^ ’ 

(34) 

2a'^R 

(35) 

G ■ 


M{R) = 

A distribution like this is expected if the bulge results from mergers. For a roughly 
constant gas fraction fg, the mass of the narrow swept-up gas shell at radius R is 


Mg{R) = 


2fga^R 


G 


SO that the shell has the equation of motion 


dt 


[Mg{R)R] + 


GMg{R)[M + M{R)] _ LEdd 
i?2 c ’ 


(36) 


(37) 


where M is the SMBH mass. From ([351 [Ml and the definition of I/Edd (equation 
[m this simplifies to 




= -2a 1 - 


M 

Ma 


(38) 


where 


M — -4 


(39) 


Multiplying through by RR and integrating once gives the first integral 


r2^2 ^ _2GMR - 2a‘^ 


1 - 


M 


R + constant 


(40) 
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For large R we have 


~ - 2 ( 7 ^ 1 


M 


(41) 


which has no solution for M < Ma- Physically this says that if the SMBH mass is 
below Ma the swept-up shell of interstellar gas cannot reach large radius because 
the Eddington thrust of the black hole wind is too small lift its weight against the 
galaxy bulge potential. The SMBH cannot remove the gas from its surroundings, 
and goes on accreting. Any gas shell it drives outwards eventually becomes too 
massive, and so tends to fall back and probably fragment. This is likely to stimulate 
star formation in the shell remnants. 

The precise value Ma depends on the average gas fraction fg. For a protogalaxy 
forming at high redshift we expect fg = 12baryon/f?matter — 0.16 (Spergel et ah, 
2003). Galaxies forming at later times may have larger fg if they have gained a lot 
of gas, or smaller fg if they have been largely swept clear of gas, or have turned a 
lot of their gas into stars. With the gas fraction fg hxed at the cosmological value 
fc = 0.16, the expression 


Ma = ~ 3.2 X lO®M0a|oo (42) 

is remarkably close to the observed relation m, even though it contains no free 
parameter. We shall see in Section 5.5 why observations tend to give an exponent 
of a slightly larger than the value 4 derived here. This agreement strongly suggests 
that SMBH growth stops at this point, although we must do some more work to 
show this (see Sections 5.2, 5.3 and 5.5 below). 

The derivation here took the simplest possible description of a galaxy spheroid 
as an isothermal sphere (cf equation I35II . We should ask if things change sig- 
nihcantly if the galaxy bulge is more complicated than this. If the potential is 
spherically symmetric but the cumulative mass M{R) is not simply linear in R, 
we still get a first integral of the equation of motion (13711 simply by multiplying 
through by M{R)R, giving the giving the condition for a swept-up momentum- 
driven shell to reach large radii. Relations very like (HI) emerge in each case, so we 
expect qualitatively similar behavior. McQuillin & McLaughlin (2012) show this 
explicitly for three widely-used density distributions (Hernquist 1990; Navarro et 
al. 1996, 1997; and Dehnen & McLaughlin 2005): the results are in practice scarcely 
distinguishable from (14211 . 

Whatever the bulge geometry, the black hole always communicates its presence 
only through the ram pressure of its wind, so we are always dealing with strongly 
radial forces in the solid angles exposed to this wind (this is not true of gas pressure, 
as we shall see in Section 5.3). It is likely that the orientation of the accretion 
disc with respect to the host galaxy changes with each new episode of accretion 
(so-called chaotic accretion. King & Pringle 2006, King et al. 2008), tending to 
isotropize the long-term effect of momentum feedback. Together with the sudden 
huge increase in the spatial scale as the critical black hole mass is reached (see the 
next Section) this may explain why the simple spherically symmetric prediction 
(14211 seems to give a surprisingly accurate estimate of the critical mass. 
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5.2 What Happens When M = Ma'? 

The result (|i^ is so close to observations that it strongly snggests that feedback 
somehow cnts off the growth of the black hole at a mass very close to this valne. 
Some feeding may continne from gas in the immediate vicinity of the hole which 
is too dense to be affected by the ram pressure of the black hole wind. A thin 
accretion disc has this property for example, but cannot have a gas mass larger 
than ~ {H/R)M ^ M withont fragmenting and forming stars. 

Bnt we still have to explain precisely how the gas is expelled. For example, 
one might worry that althongh momentnm—driving can push the ISM away and 
inhibit central accretion, some kind of infall and SMBH accretion might restart 
shortly after momentnm-driving is switched off, perhaps leading to alternating 
stages of quiescence and growth, eventually to masses far above Ma- Observations 
show that black hole accretion occurs preferentially in gas-rich galaxies (cf Vito 
et al. 2014), so it seems that the black hole must largely clear the galaxy bulge of 
gas to terminate its growth. We will see later (Section 7) that if no other process 
than momentum—driving operated, this requirement would indeed lead to black 
hole masses significantly larger than Ma, in conflict with observation. 

This last point means that the way black hole growth influences the host galaxy 
must change radically when M = Ma- It is straightforward to see why it should. We 
saw above that for M < Ma the Eddington thrust cannot push the wind shocks to 
large R. As a result the wind shock remains efficiently Compton cooled, enforcing 
momentum-driving. It follows that the SMBH cannot stifle its own growth if 
M < Ma. But all this changes once the SMBH exceeds the critical mass (I42II . 
Now even for a very small increment {0{Ri^i/Rc) ~ 10“^) of M above Ma, a 
momentum-driven shell can reach the critical radius Re- Crucially, this means 
that the wind shocks are no longer efficiently cooled: they become energy-driven. 
The shocked wind gas can now use all of its energy to push the interstellar gas as 
it expands into the host bulge. This motion becomes explosive and rapidly reaches 
kiloparsec lengthscales, comparable with the size of the bulge itself, rather the 
much smaller (parsec) scales of the momentum-driven phase. 

So the real significance of the M — a relation is that it marks the point where 
outflows undergo a global transition from momentum- to energy—driving. 


5.3 Clearing Out a Galaxy: the Energy-Driven Phase 

We know from (O that an energy-driven outflow has more than enough energy 
to remove the interstellar gas entirely, and so presumably suppress further SMBH 
growth. Here we examine how this works in detail. 

Once M > Ma the outflow geometry changes completely (see Eigure 8). The 
shocked wind region is no longer narrow (as in the upper panel of Figure 8), but 
large and expanding because of its strong thermal pressure (lower panel of Figure 
8). The shocked wind’s thermal expansion pushes its shock inwards where it must 
hover at the cooling radius Rq (Zubovas et al. 2013). If it tries to move within 
Rq, momentum driving instantly pushes it out again (remember M > Ma). 

The shocked wind rapidly evens out its internal pressure as it expands at its 
sound speed ~ 0.03c, so we take this pressure P as uniform over this region (but 
changing with time). The contact discontinuity at the outer edge of the shocked 
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wind sweeps up the surrounding shocked ISM as before, but now has the equation 
of motion 


_d 

dt 


Mg{R)R 




( 43 ) 


where the pressure P is much larger than the ram pressure appearing in (|371l. 
In the second term on the Ihs we have neglected the contribution GMgM/R^ of the 
black hole gravity, as R^ Rq > -Rinf • To make the problem determinate we need 
the energy equation. (This did not appear explicitly in the momentum-driven case 
because it was equivalent to the defining condition that all the wind energy not 
associated with the ram pressure was rapidly lost to radiation.) Here the energy 
equation constrains the pressure P by specifying the rate that energy is fed into 
the shocked gas, minus the rate of PdV working on the ambient gas and against 
gravity: 


d 

'AttR^ 3„ 

- 


'47r 3 ' 

dt 

3 '2^ 

2 -^Edd 

dt 


-4/s- 


(44) 


We take a specihc heat ratio 7 = 5/3, use (1231) for the energy input from the 
outflow and (1351) to simplify the gravity term GM{R)M{R)/R?. Now we use (1431) to 
eliminate P from (1441) , and replace the gravity terms as before using the isothermal 
expression for M{R). We take the AGN luminosity as IL^dd to allow for small 
deviations from the Eddington value. This gives 




M{R)R 




M{R)R 


+ 2/g-^R 


(45) 


and so 

I^TEdd = + 3RM+ + Wfg^R. (46) 

This describes the motion of the interface (‘contact discontinuity’ in the lower 
panel of Figure 8 ) between wind and interstellar gas in the energy-driven case, 
replacing equation dSIl) in the momentum-driven case. 

The energy-driven regime applies as soon as the SMBH mass reaches Ma, and 
we will see that the host ISM is now quickly removed. We assume M = Ma in 
■bEddj and see that (1461) has a solution R = Vet with 


2r]lc = 3^ + 10«e 


(47) 


The assumption Ve << u leads to a contradiction {ve ^ 0.01c >> a), so 


Ve ~ 


2r]la^c 


1/3 


~ 925;^/^(J2^o km 


(48) 


This solution is an attractor. Figure 9 shows that all solutions quickly converge 
to it, regardless of initial conditions. Physically, its meaning is that if shock cooling 
is ineffective, the extra gas pressure accelerates the previously momentum-driven 
gas shell to this new higher velocity. Figure 9 also confirms that if the driving by 
the AGN switches off when the contact discontinuity is at radius Rq, it decelerates 
as predicted by the analytic solution of (1461) with I/Edd = 0 found by King et al. 
( 2011 ): 









( 49 ) 
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where x = R/Rq > 1. Noting that Ve depends only weakly (as Ve ~ on 

the luminosity, these results show that fluctuations - or even the intermittent 
disappearance - of the AGN luminosity have almost no effect on the outflow once 
it has started, because the flow still has a large reservoir of thermal energy available 
for driving. In particular an outflow can persist long after the central AGN has 
turned off, and the real agency driving an observed outflow may have been an 
AGN even if this is currently observed to be weak or entirely absent. 

The solutions (HSll or describe the motion of the contact discontinuity 
where the shocked wind encounters swept-up interstellar gas (see Figures 8 and 
9). This interface is strongly Rayleigh-Taylor (RT) unstable, because the shocked 
wind gas is highly expanded and has much lower mass density than the swept-up 
interstellar gas outside it, so that we have a light fluid underneath a heavy one. 
The RT instability leads to strong overturning motions even on small scales, and 
so is difficult to handle numerically. Deductions concerning the mean velocity and 
energy of the outflow, and its average spatial scale R{t), are likely to be believable, 
and agree closely with observations (see below) but we should be very cautious 
about results depending strongly on the detailed nature of the interface between 
the shocked wind and the swept-up interstellar medium. The RT instability is 
probably the reason that the high-speed (~ 1000 kms“^) outflows with prodi¬ 
gious mass rates we predict here are generally seen with much of the outflowing 
gas in molecular form. Apparently the interstellar gas entering the forward shock is 
efficiently cooled by two-body radiation processes. A preliminary analysis (Zubo¬ 
vas & King 2014) suggests that the interstellar gas overtaken by the forward shock 
is likely to have a multiphase structure. Most of it cools all the way from the shock 
temperature ~ 10^ K back to low temperatures, ending in largely molecular form, 
even though it is entrained in an outflow with the ~ 1000 kms“^ velocity of the 
forward shock. But cooling is affected by the topology of the gas flow and the total 
area of interfaces between different gas phases. A full numerical calculation of this 
is currently impossible, so for the time being we can only make comparison with 
simple estimates, as here. 

The mass outflow rate is fixed by how fast the outer shock overtakes the ISM 
and entrains new interstellar gas ahead of the contact discontinuity. The ISM 
ahead of the shock is at rest, so this runs into it at a speed giving a velocity jump 
by a factor (7 -|- l )/(7 — 1 ) in the shock frame (where 7 is the specific heat ratio: 
see e.g. Dyson & Williams 1997 for a derivation). This fixes its velocity as 

vout = ^ 1230 ( 72/0 ' km s-i (50) 

(where we have used 7 = 5/3 in the last form, and /c 0.16 is the cosmological 
value of /g). This implies a shock temperature of order 10^ K for the forward (ISM) 
shock (as opposed to ~ 10^°“^^ K for the wind shock). Since the outer shock and 
the contact discontinuity were very close together as energy-driven flow took over 
from momentum-driven flow (see Figure 8 ) this means that the outer shock is at 


Rout(t) = '^^4r^R{t) = 


This gives the mass outflow rate as 


. _ dM(Rout) _ (7+l)/g'^% 

- -d^-- -G- 


(51) 


G 


( 52 ) 
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Radius / kpc 

Fig. 9 Evolution of an energy—driven shock pattern for the case u = 200 kms~^, fg = 10“^ 
computed numerically from the full equation gg. Top: radius vs time, middle: velocity vs 
time, bottom: velocity vs radius. The curves refer to different initial conditions: black solid - 
Rq = 10 pc, no = 400 kms“^; blue dashed - Rq = 100 pc, vq = 1000 kms~^; red dot—dashed 
- Rq = 50 pc, uq = 200 kms“^. All these solutions converge to the attractor gg. The vertical 
dashed line marks the time t = 10® yr when (for this case) the quasar driving is switched off. 
All solutions then follow the analytic solution gg. A case where the quasar remains on for a 
Salpeter times ~ 4 x 10^ yr would sweep the galaxy clear of gas. (From King et ah, 2011) 
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For comparison the mass rate of the black hole wind, assuming M = Ma, is 


Mw = mMEdd — 


4/cmg^ 

rjcG 


(53) 


This is much smaller than the outflow rate Mout it drives, so we define a mass¬ 
loading factor as the ratio of the mass flow rate in the shocked ISM to that in the 
wind: 


/l 


-^out 

Mw 


^(7 + 1 ) fg Rc 
4m fc (J^ 


(54) 


Then we have 




vh + 1 ) fg Rc 
4 /c<t2 


MEdd- 


(55) 


If the AGN radiates at a luminosity ~ Z/Edd, we have R = Ve, and (1481) gives 


/l 



fg 


2/3 ^ 1/3 


~ 460cr 


- 2/3 

200 


m ’ 


(56) 


and 

Mout 4060 ( 72 ^°/^/^/^ M 0 yr-i (57) 

for typical parameters, /g = fc and 7 = 5/3. The total gas mass in the bulge is 
roughly Mg ~ lO^/gMo- (from equation [T|) . Clearly if the outflow persists for a 
time tciear ~ Mg/Mont ~ 1 X yr it will sweep away a large fraction 

of the galaxy’s gas. The precise outflow duration needed for this depends on both 
the type and the environment of the galaxy, in practice leading to three parallel 
but slightly offset M — a relations (see Section 5.5 below). 

Equations (ISOllFfl) give 


iMwu" ~ ^Montvlof (58) 

So most of the wind kinetic energy ultimately goes into the mechanical energy of 
the outflow, as we would expect for energy driving. The continuity relations across 
the contact discontinuity show that if the quasar is still active, the shocked wind 
retains 1/3 of the total incident wind kinetic energy M^v^ j2, giving 2/3 to the 
swept-up gas represented by Mout- 

Equation (1581) means that the swept-up gas must have a scalar momentum 
rate greater than the Eddington value 7/Edd/c, since we can rewrite it as 


2Mw 2Mout 

where the P’s are the momentum fluxes. With Pw = TEdd/c, we have 

• 1 /2 

D D /Mout\ T/Edd j.1/2 on^Edd -1/3,1/6 

Pout = Pw ^— = - /t 20 -(Tjnn ‘ • 

V Mw / c C 200 


(59) 


(60) 


Observations of molecular outflows consistently show MoutVout > 7/Edd/c, and in 
particular Cicone et al. (2014) find that momentum rates 20P/c are common. This 
is an inevitable consequence of mass-loading (/^ > 1). These high momentum rates 
are important, as they are probably the way that the galaxy resists the accretion 












AGN Outlflows and Feedback 


31 


that cosmological simulations suggest still continues at large scales (Costa et al. 
2014). 

Recent infrared observations show abundant evidence for molecular outflows 
with speeds and mass rates similar to (1501) and (1571) . Feruglio et al. (2010), Rupke 
& Veilleux (2011) and Sturm et al. (2011) find large-scale (kpc) flows with Vout ~ 
1000 kms“^) and Mout ~ 1000 MQyr“^) in the nearby quasar Mrk 231. Other 
galaxies show similar phenomena (cf Lonsdale et al. 2003, Tacconi et al. 2002, 
Veilleux et al. 2009 Riffel & Storchi-Bergmann (2011a, b) and Sturm et al. 2011: 
see Tables 1 and 2 of Zubovas & King 2012a for a detailed comparison with the 
theoretical predictions). In each case it appears that AGN feedback is the driving 
agency. There is general agreement for Mrk231 for example that the mass outflow 
rate Mout and the kinetic energy rate Eout = MoutVout/S are too large to be driven 
by star formation, but comparable with values predicted for AGN feedback. 

It appears that energy-driven outflows from SMBH which have just reached 
their M — a masses should be able to sweep galaxy spheroids clear of gas. A robust 
observational test of this is the expected mechanical luminosity (cf equation 1231) 

Tmech ~ I^Edd ^ 0.05L, (61) 

where L = IL^dd is the observed AGN luminosity. This is investigated by Gicone et 
al. (2014). As their Figure 12 shows, observation does largely confirm the relation 
EB). If the AGN are close to their Eddington luminosities (so that L ex M ex 
and I ~ 1), the clearout rate oc feauation l57l) should scale linearly with the 

driving luminosity L. Figure 9 of Cicone et al. shows evidence for this correlation, 
with normalization close to that predicted. 


5.4 Effects of a Galaxy Disc: Stimulated Star Formation and Outflow Morphology 


We have so far discussed galaxy spheroids in isolation. This is in line with the 
observational evidence (see Kormendy & Ho 2013 for a review) that the SMBH 
scaling relations apply only to this component of a galaxy, and are essentially 
unaffected by the presence of a galaxy disc. In particular we suggest that the 
critical M — a black hole mass is set by small-scale momentum-driven outflows 
interacting with only a very small central part of the bulge. But the energy-driven 
outflows we considered in the last subsection are global: they expand to far greater 
scales, and unless the galaxy is an elliptical must inevitably encounter its disc as 
they expand. In a gas-rich galaxy the gas in the innermost disc at radius Rq must 
be close to self-gravitating. Assuming that the potential is roughly isothermal, 
it is straightforward to show that this implies a gas density ~ Pd ~ /Rq, i.e. 
greater than the bulge gas density by the factor ~ l//g ~ 10. We see from equation 
(|i6l) that higher gas densities mean lower spherical outflow velocities, as they meet 
greater resistance. So when an initially spherical outflow encounters a high-density 
gas disc it flows around it, over its plane upper and lower faces. But the pressure 
in the outflow is at least initially far higher than in the disc: we can read off the 
pressure at the contact discontinuity from equation (|i3l) as 


-PcD 


fga\2e7^ + v^e) 
leGR^ 


r 2 2 

7rGR2 


(62) 





32 


King and Pounds 


and estimate the pressure at the forward shock into the ISM as 

D 2 2/ga Ve 2 i'ao\ 

P,s = ^p{R)ve = C. gPcD. (63) 

By contrast the mid-plane pressure in a disc close to self-gravitating is 

4 

-Pdisc ~ pci ~ ~ 

where we have assumed the sound speed Cs ~ cr and the self-gravity condition 
Gp ~ with Q = \f2alR^ the Kepler frequency at disc radius R^. Thus when 
the outflow shock arrives sX R = R^ its pressure is a factor ~ (yejo')^ ~ 25 larger 
than the disc’s, and this remains true until the outflow shock has travelled out 
to radii R > R^Ve/u ~ ^Rd- Any such compression must trigger a burst of star 
formation in the disc (cf Thompson et al. 2005, Appendix B), and here it rises to 
values 

i;* ~ 20006_3<J2^S^^/^/^i?kpe M0 kpc-2 (65) 

(Zubovas et al. 2013), where e* = 10“^e_3 is the efficiency of massive stars in 
converting mass into radiation, and we have substituted for Ve using (I48II . Zubovas 
et al. (2013) show that this leads to a starburst of total luminosity 

L* ~ 5 X ergs"\ (66) 

where Z /46 is the AGN luminosity in units of 10^® ergs“^. Such systems would 
appear as ULIRGs. 

This suggests that in a galaxy with both a bulge and a disc, the clearout phase 
leaves the galaxy bulge without gas, but may be accompanied by a starburst in 
the disc. Recent observations of dusty QSOs appear to show this, with the black 
hole mass already on the M — relation o, and so fully grown (Bongiorno et 
al. 2014). In an elliptical on the other hand, clearout must leave the galaxy ‘red 
and dead’. 

Since a galaxy disc is a major obstacle to an outflow, it follows that it may be 
able to divert a quasi-spherical outflow into a bipolar shape. This is particularly 
true in cases where the SMBH mass grows only a little, in a minor accretion 
event. Zubovas et al. (2011) suggest that the gamma-ray emitting bubbles disposed 
symmetrically about the plane of the Milky Way (Su et al. 2010) may be the 
remnants of a relatively recent and rather weak event like this. 


5.5 The Three M — a relations 

So far in this Section we have seen that the arrival of the black hole mass at the 
M — a relation means that its feedback makes a radical change from momentum¬ 
driving to energy-driving. The energy-driven phase which clears the gas from a 
galaxy bulge is short and violent. But it is clear that the black hole must inject 
a non-negligible amount of energy to eject the gas, and this requires accretion 
energy, i.e. some black hole mass growth. Evidently if the mass increment AM 
needed for this is S> Mcr we will have failed to explain the M — a relation. 

The mass increment AM is influenced by two factors. First, it must require 
significantly less SMBH growth to remove the gas from a spiral galaxy with a 
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Velocity dispersion, km/s 


Fig. 10 The four (in reality three, as cluster spirals are rare) M — a relations (solid lines) and 
their combined effect on observational fits (dashed line). All solid lines have slopes M oc 
and the dashed line has M cc . The grey area is the approximate locus of data points in 
Figure 3 of McConnell et al. (2011). (From Zubovas & King, 2012b) 


relatively small bulge, than for example an elliptical, where the much larger bulge 
mass means that energy-driving by the central SMBH wind must continue for 
longer in order to expel the remaining gas. Zubovas & King (2012b) find that 
energy-driving, and therefore SMBH mass growth above M^, must continue only 
for about 4 Myr (about 0.1 Salpeter times) in a typical spiral, but for for about 2 
Salpeter times in an elliptical. So the final SMBH mass in a spiral is close to Ma 
but in an elliptical it can reach 


Mfinai ~ ~ 7.5M,, (67) 

The second factor affecting M,Mi, is the galaxy environment. Cluster ellipticals 
can gain gas as they orbit through the intracluster gas. Some Brightest Cluster 
Galaxies (BCGs), which are near the centre of the cluster potential, are known to 
contain unusually massive SMBH (McConnell et al. 2011). Taking account of the 
extra black hole mass growth required to remove the bulge, and the mass a galaxy 
may gain from its surroundings, implies three parallel but slightly offset M — a 
relations for spirals, field and cluster ellipticals (see Fig 10). In principle there is 
also a relation for cluster spirals, but these are rare. We see from the Figure that 
the spread in offsets means that an observed sample drawn from galaxies of all 
three types would tend to produce a slope slightly bigger than the individual ones 
for each type, perhaps accounting for the slight discrepancy between the observed 
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Fig. 11 The four (in reality three, as cluster spirals are rare) M — a relations (solid lines) and 
their combined effect on observational fits (dashed line). All solid lines have slopes M oc 
and the dashed line has M oc cr^. The grey area is the approximate locus of data points in 
Figure 3 of McConnell et al. (2011). (From Zubovas & King, 2012b) 


overall slope a = 4.4±0.3 and the theoretical value of 4. All three types of galaxies 
obey a similar M — relation © within the errors, as growth of the SMBH above 
Ma goes together with higher 


6 THE SMBH - BULGE MASS RELATION 

6.1 Feedback and the M — relation 

In the Introduction we noted the observed proportionality o between M and 
Ml, as well as the M — a relation. So far we have concentrated almost entirely on 
the second of these relations, and suggested that it arises because the black hole 
feedback itself directly limits the mass reservoir available for black-hole growth. 
Quite independently of details, almost every discussion of this relation adopts this 
view (see Section 7 below). 

But the character of the M — Mi, relation must be very different. Since we are 
assuming that feedback ensures that the black hole mass M is set by a we cannot 
argue that M is independently set by Mi,. But reversing the argument to suggest 
that the black hole mass M sets Mi, is also implausible, since is in the form of 
stars. 
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So there can be no directly causal connection between the black hole mass M 
and the stellar bulge mass (Indeed one view - see Section 6.3 below - asserts 
that the connection is purely statistical.) Instead, their relation must arise because 
whatever determines makes it proportional to cr^. Empirically, we already know 
that this is approximately true, at least for elliptical galaxies, the largest spheroids 
of all, because these are observed to obey the Faber-Jackson (1976) relation 

L, ~ 2 X lO^°L0cr|oo, (68) 


Here L* is the total stellar luminosity and mass of an elliptical, so for mass-to-light 
ratios ~ 5 we immediately get the stellar mass as 

M* ~ 1 X IO^^MqctIoo ~ lO^M (69) 


There is now general agreement that this relation, like the M — a relation, may re¬ 
sult from feedback inhibiting and ultimately suppressing the process that produces 
it. The difference is that here the feedback is from stars, and what ultimately has 
to be suppressed is star formation. Several papers make this point, starting with 
Murray et al. (2005). Power et al. (2011) show that this approach gives a bulge 
stellar mass 


Mb 


O.UfgtHa^ 
6 *; CG 


(70) 


where e* ~ 2 x 10~^ measures the total luminous energy yield from a main- 
sequence star in terms of its rest-mass energy M*c^, and tu is the Hubble time. 
Comparing with (I39II we get 


~ W-^Mb, (71) 

which is similar to observational estimates (cf equation [1} . Both the M — a and 
the M — Mb relations hold for elliptical galaxies, so equation (EOl) automatically 
reproduces the Faber-Jackson (1976) relation for typical mass-to-light ratios. In 
this view, the similarity of the SMBH and stellar (Faber-Jackson) M,Mb oc. 
relations (HHMI follows directly because both result from momentum feedback, 
and the ratio M/Mb ~ 10“^ reflects the relative efficiencies of the black hole and 
stellar versions. 


6.2 The M — a relation for Nucleated Galaxies 

A similar argument (McLaughlin et al. 2006) shows that for nucleated galaxies 
(i.e those whose central regions are dominated by nuclear star clusters, with no 
detectable sign of the presence of a supermassive black hole) there should be an 
offset M — a relation between the mass of the cluster and the velocity dispersion, 

i.e. 

Me ~ 20Ma ~ 6 X loVlooM©. (72) 

Typically these galaxies are small, with cr < 120 kms“^. The factor ~ 20 offset 
in cluster mass for a given cr arises because momentum-driving by an ensemble of 
cluster stars is about 20 times less efficient per unit mass than from a black hole 
accreting at the Eddington rate. 
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6.3 Mergers and the M — Mf, relation 

Jahnke & Maccio (2011) offer a radically different interpretation of the M — 
relation. Building on earlier work of Peng (2007) they assume that black hole and 
bulge masses are built up by repeated mergers of smaller galaxies with uncorrelated 
M and Mf,. They follow this evolution using dark matter halo merger trees, and 
as a result of the central limit theorem find that M is roughly proportional to M^, 
with the scatter decreasing for larger masses, where there have been more mergers. 
They conclude that the SMBH - bulge scaling relation may have an explanation 
that is largely or even entirely non-causal. 

But it is hard to accept that there is no more physics in the SMBH scaling 
relations than this. First, the actual ratio M/Mj, is left undetermined by this pro¬ 
cedure. Second, to get from the M — M^ relation to M — a requires one to assume 
something like the Mf, oc relation (IMl) implied by Faber-Jackson, so physics 
presumably must enter here too (cf the subsection above). Finally, it would seem 
a remarkable coincidence that the outcome of this indirect process by chance pro¬ 
duces an M — a relation exactly equivalent to requiring that the SMBH Eddington 
thrust should just balance the weight of the bulge gas. 


7 MOMENTUM, ENERGY OR RADIATION? 

The study of AGN outflows and their effects on the host galaxy has two main 
aims. A viable picture must explain both the scaling relations, and simultaneously 
the fact that galaxy spheroids appear ultimately to be largely swept clear of gas 
by high-speed molecular outflows which have significantly greater scalar momenta 
Mv than the radiative value Ljc of the central AGN (the clearout problem). The 
discussion given above offers plausible physical grounds that the shock interaction 
characterising the black hole wind feedback changes from momentum-driven, act¬ 
ing on small spatial scales near the black hole, to energy-driven, instead acting 
globally on the whole galaxy bulge and producing a high-energy clearout of its 
gas. The M — a relation marks the point where this transition occurs in a given 
galaxy. We will argue below (Section 8) that observations support this picture of 
local-global transition in several ways, but before accepting this conclusion we 
should consider other possibilities. 

First, the switch from momentum to energy-driving depends on the details of 
gas cooling. It is sometimes argued (e.g. Silk & Nusser 2010) that strong cooling 
of the ambient interstellar medium enforces momentum—driving by the central 
SMBH throughout. In fact cooling the ambient gas is not relevant to the question 
of energy or momentum driving: as we have seen, it is the cooling of the shocked 
wind gas which decides this. But as we emphasised in Section 4, at least some 
of the physics of the suggested momentum-energy switch is still beyond a full 
numerical treatment with realistic parameters. It is sensible then to check our 
treatment above by considering the momentum-driven and energy—driven cases in 
isolation, and then the effect of direct radiation pressure. 
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7.1 Wind Momentum Driving 


We first simply assume that a black hole wind acts on its surroundings by pure 
momentum—driving alone, at all radii. For generality we let the pre-shock wind 
have speed Vw and take its mechanical luminosity MwVw/2 as a fixed fraction a of 
-f'Edd: i-G. we do not explicitly assume that the wind has the Eddington momentum, 
as seems to hold for UFOs. Then the momentum feedback first becomes important 
at a critical black hole mass Merit roughly given by equating the wind thrust 
MwVw = 2aL^dd/v fo the weight 


_ GM{R)Mg{R) 


G 


(73) 


of the overlying gas in an isothermal potential (cf equation I37II . With I/Edd = 
4GMcritc/«: we get 


Merit 


Vw 

2ac 


Ma. 


By definition a = Mwv'^/2L-e^^ and I/Edd = so 


Merit 
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(75) 


(cf Fabian 1999). We see that for general wind parameters the critical mass differs 
from Mcr. We find Merit = Ma only if Vw = ycMEdd/Mw, which is equation (1^ . 
This immediately implies that the wind momentum is Eddington, i.e. MwVw = 
yM^ddC = 7/Edd/c. In other words, assuming pure momentum driving gives the 
critical mass as Ma if and only if the driving wind has the Eddington momentum, 
i.e. has the properties observed for UFOs. 

But pure momentum driving is unable to drive off the bulge gas without a sig¬ 
nificant increase of the black hole mass above Merit- Several authors have reached 
this conclusion (cf Silk & Nusser 2010, McQuillin & McLaughlin 2012). More¬ 
over, if galaxy bulges accrete at the rates suggested by cosmological simulations 
it seems unlikely that any hypothetical momentum-driven outflows would have 
enough thrust to prevent infall and so could not turn off star formation defini¬ 
tively (cf Costa et al. 2014). We conclude that pure momentum-driving, even 
given the lack of a likely shock cooling process, probably does not give a realistic 
picture of the interaction between SMBH and their hosts. 


7.2 Wind Energy Driving 

The direct opposite case from that considered in the last subsection is pure energy¬ 
driving by winds, where radiative cooling is assumed negligible throughout. This 
was often the implicit assumption in early treatments (e.g. Silk & Rees 1998, 
Haehnelt et al. 1998). The equation of motion for this case is (1461) . This shows 
that gas is always driven out at constant speed for any SMBH mass, however 
small: setting R = Vet and using the definition of I/Edd (cf equation 1 171) gives the 
speed Ve as 

3 nG^crjM 
Sfgua'^ 


(76) 
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This expresses the fact that the adiabatically expanding shocked gas pushes the 
interstellar gas away as a hot atmosphere for any SMBH mass. One can easily find 
the corresponding mass outflow rate by setting Moutw|/2 ~ Lw, since we know 
that the outflow mechanical luminosity is a significant fraction of the driving wind 
mechanical luminosity Lw = rjlL^^^/2. 

To dehne a critical SMBH mass for energy-driven outflow one has to impose 
a further condition. This is usually taken as We ~ u, defining some kind of escape 
velocity. But it is not obvious that this is appropriate: the outflow is driven by 
pressure, so the ballistic escape velocity is not relevant. Even if the AGN switches 
off, the residual gas pressure still drives outflow for a long time (cf Fig. 9). If we 
nevertheless impose this condition we find a critical mass 

Menergy = = —Me = 0.02M^ = 6 X lO®M0Crfoo (77) 

TrG^jyc 77 c 

which is a factor “iaI'qc ~ 1/50 too small in comparison with observations. 

Silk & Rees (1998) considered the growth of a protogalaxy (i.e. gas with fg ~ 1) 
around a supermassive seed black hole which formed earlier, but their argument 
applies to the coevolution of the SMBH and host also, provided we take fg ~ 0.1. 
They assume the wind sweeps mass into a shell with speed Vs, and implicitly 
neglect pressure work, and the fact that energy is shared between the shocked 
wind and the swept—up outflow. This would imply a relation 

Lw = 2-Kr‘^ fgp{r)vl, (78) 

as each new shell of mass 47rr^p(r)i!s now simply acquires kinetic energy vl/2 as 
it is swept up. Using the isothermal relation pill and requiring Vs ~ u gives 

where fw = Tui/T/Edd- The neglect of pressure work overestimates the wind-driving 
efficiency, so this mass is even smaller than (EZI). It is clear that wind energy¬ 
driving of this type does not correctly reproduce the M — a relation, giving a 
critical mass too low by factors 50 - 100. 

A more promising approach has recently been outlined by Nayakshin (2014), 
Zubovas & Nayakshin (2014) and Bourne et al. (2014), who consider the effects of 
strong inhomogeneity of the bulge gas. They assume first that inverse Compton 
shock cooling may not be effective because of two-temperature effects (but see 
Section 4.2 above). Second, they suggest that sufficiently dense clouds of interstel¬ 
lar gas would feel a net outward force ~ pv^ per unit area when overtaken by a 
free-streaming UFO wind of preshock density p and speed v, thus mimicking a 
momentum-driven case. The density of these clouds is a factor 1/fg ~ 6 below 
the star-formation threshold. If most of the ISM gas mass is in the form of such 
clouds, SMBH feeding must stop when the outward force overcomes gravity, which 
gives a relation like (13911 up to some numerical factor. This idea throws up several 
gas-dynamical problems. First, a cogent treatment must explain the origin and 
survival of clouds at densities close to but just below the star formation threshold, 
which must contain most of the interstellar gas. The clouds must be completely 
immersed in the wind, so the net outward force on them is a surface drag, which 
is dimensionally also ~ pv^ per unit area. Estimating this surface drag requires 
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knowledge of how the cloud-wind interfaces evolve on very small scales. Since 
these are formidable theoretical tasks, we should ask for observational tests. The 
main difference from the quasispherical momentum—driven case is that instead 
of being radiated away, most of the energy of the UFO wind now continuously 
drives the tenuous intercloud part of the ISM out of the galaxy at high speed. 
If this tenuous gas is a fraction ft of the total gas content, equations (1501 EH) 
show that for SMBH masses not too far from Ma this outflow should have speed 
«out ~ kms“^ and mass-loss rate Mout ~ 4000/t M 0 yr“^, and so be 

potentially observable for many AGN spheroids. From the work of Section 5.4 one 
might also expect a continuously elevated star formation rate in the central parts 
of their galaxy discs also, which is not in general observed. 


7.3 Cosmological Simulations 

Cosmological simulations often produce an empirical M—a relation as part of much 
larger structure formation calculations. Limits on numerical resolution inevitably 
require a much more broad-brush approach then adopted here. The effect of the 
SMBH on its surroundings is usually modelled by distributing energy over the gas 
of the numerical ‘galaxy’ at a certain rate. This injected mechanical luminosity 
is then iterated until the right relation appears. This empirical approach (e.g. di 
Matteo et al. 2005) seems always to require a mechanical luminosity 0.051/Bdd to 
produce the observed M — a relation. This is precisely what we expect (cf equation 
1231) for a black hole wind with the Eddington momentum Moutv = I/Edd/c. 

But the success of this procedure is puzzling. If the ambient gas absorbed the 
full numerically injected mechanical luminosity 0.051/Edd the resulting outflows 
would give the energy-driven (13211 or Silk-Rees mass (17911 above, which are too 
small compared with observations. The fact that cosmological simulations instead 
actually iterate roughly to the observed M — a value dH suggests that they some¬ 
how arrange that the injected energy only couples to the gas at the very inefficient 
rate which occurs in momentum driving, or possibly that the numerical gas dis¬ 
tribution is highly inhomogenous. The real physics determining this in both cases 
operates at lengthscales far below the resolution of any conceivable cosmological 
simulation, so the inefficiency must be implicit in some of the ‘sub-grid’ physics 
which all such simulations have to assume (cf Costa et al. 2014, Appendix B). 


7.4 Radiation Driving 

l.f.l Electron scattering opacity 

We remarked in the Introduction that in principle direct radiation pressure is 
the strongest perturbation that a black hole makes on its surroundings, but its 
effects are more limited than this suggests. As we already suggested, the reason is 
that in many situations radiation decouples from matter before it has transferred 
significant energy or momentum. This is particularly likely for radiation emitted 
by an AGN in the center of a galaxy bulge. The gas density (cf equation 1341) is 
sharply peaked towards the center, and so is its tendency to absorb or scatter the 
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radiation from the accreting black hole. The electron scattering optical depth from 
gas outside a radius R for example is 

r{R) = l^ Kp{r)dr=^. (80) 

which is mostly concentrated near the inner radius R. This means that gas initially 
close to the AGN is probably swept into a thin shell by its radiation, and so at 
radius R has optical depth 

“ nGR ’ 

very similar to the undisturbed gas outside R (cf equation 1801) . Gas distributed 
in this way has large optical depth near the black hole when its inner edge R 
is small (i.e. less than the value Rn specihed in equation [82] below). Then the 
accumulating accretion luminosity L of the AGN is initially largely trapped and 
isotropized by electron scattering, producing a blackbody radiation held whose 
pressure grows as the central AGN radiates. This growing pressure pushes against 
the weight W = feauation l73ll of the swept-up gas shell at radius R. This 

is exactly like the material energy-driving we discussed in Section 5.3, except here 
the photon ‘gas’ has 7 = 4/3 rather than 7 = 5/3 there. Clearly the effectiveness 
of this radiation driving is eventually limited because the shell’s optical depth falls 
off like 1/7? as it expands. The force exerted by the radiation drops as it begins 
to leak out of the cavity, until for some value rtot(7?) ~ 1 it cannot drive the shell 
further. 

This shows that the sweeping up of gas by radiation pressure must stop at a 
‘transparency radius’ 

(82) 

where (up to a logarithmic factor) the optical depth rtot is of order 1 , so that the 
radiation just escapes, acting as a safety valve for the otherwise growing radiation 
pressure. This process is discussed in in detail by King & Pounds (2014), who 
suggest that the stalled gas at 7?tr may be the origin of the ‘warm absorber’ 
phenomenon (cf Tombesi et al. 2013). The radius Rtr is so small that very little 
accretion energy is needed to blow interstellar gas to establish this structure, and 
to adjust it as the galaxy grows and changes a. 

7 . 4.2 Dust 

At larger radii much of the cold diffuse matter in the galaxy bulge may be in the 
form of dust. The absorption coefficient of dust depends strongly on wavelength 
and is far higher than electron scattering in the ultraviolet, but decreases sharply 
in the infrared (e.g. Draine & Lee, 1984). The energy of an ultraviolet photon 
absorbed by a dust grain may be re-emitted almost isotropically as many infrared 
photons, which then escape freely. The net effect is that dusty gas feels only the 
initial momentum of the incident UV photon, while most of the incident energy 
escapes. Then a spherical shell around an AGN would experience a radial force 
~ L/c, where L is the ultraviolet luminosity, as long as it remained optically 
thick to this kind of dust absorption. This is dynamically similar to wind-powered 
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flows in the momentum-driven limit, and this type of radiation—powered flow is 
often also called ‘momentum-driven’, even though the physical mechauism is very 
different. 

An important distinction between the wind and radiation-powered cases is that 
ambient gas in the path of a wind cannot avoid feeling its effects, whereas this is 
not true for radiation, as the gas may be optically thin. Galaxies are generally 
optically thin to photons in various wavelength ranges, and a radiation-driven 
shell may stall at finite radius because its optical depth r becomes so small that 
the radiation force decouples, as we saw in the electron-scattering case. Ishibashi 
& Fabian (2012, 2013, 2014) appeal to this property to suggest that star formation 
in massive galaxies proceeds from inside to outside as radiation-momentum driven 
shells of dusty gas are driven out and then stall at the dust transparency radius 
fl^dust — («d/'^)^tr. For large dust opacities ~ IOOOk this can give i?dust ~ 
50 kpc. In contrast galaxies are probably never ‘optically thin’ to winds, and the 
density of a black hole wind is always diluted as 1/i?^, so it inevitably shocks 
against a swept-up shell of interstellar gas at large R. 

The mathematical similarity (cf eg 1371) between wind-powered and radiation- 
powered momentum driving allows an empirical estimate of an M — cr relation 
for the latter if we assume that observed AGN define the relation, and that 
their observed luminosities correspond to L/L^^dd ~ 0.1 — 1. This gives Merit = 
{L-^dd/L)Mcr ~ 1 — IOMct (Murray et al. 2005). Optical depth effects might narrow 
this range closer to the observed one (Debuhr et al. 2012). This suggests that radia¬ 
tion driving might be compatible with the M — a relation, but a momentum-driven 
outflow like this can never simultaneously reproduce the high-velocity molec¬ 
ular outflows characterising the clearout phase. In particular its momentum is 
Ljc < I/Edd/c, considerably smaller than the observed ~ 20LEdd/c of such flows 
(see Section 5.3). In other words, we have the usual difficulty that momentum¬ 
driving can accommodate the M — a relation, but not simultaneously solve the 
clearout problem. 

One way of possibly overcoming this (e.g. Faucher-Giguere et al. 2012; Faucher- 
Giguere & Quataert, 2012) is to assume (cf Roth et al. 2012) that instead of de¬ 
grading incident high-energy photons to lower-energy ones that escape freely, the 
effect of dust absorption is to retain much of the incident radiant energy. Then if 
the dust is distributed spherically and is in a steady state the radiation force on it 
is tL/c, where r is the radial optical depth of the dust (cf Roth et al. 2012). This 
form of radiation driving of optically thick dust can in principle produce outflows 
whose scalar momenta are boosted above that of the driving luminosity L/c by 
a factor ~ r because photons may be reabsorbed several times. For r ^ 1 the 
radiation field is effectively trapped and presumably approaches a blackbody form 
(cf the discussion of the electron scattering case above), limiting the boost. 

Evidently for radiation driving of dust to solve simultaneously both the SMBH 
scaling and clearout problems requires a sharp transition in the properties of the 
dust opacity at the critical M — a mass. This must change the outflows from effec¬ 
tively momentum-driven (incident photons are absorbed but their energy escapes 
as softer photons) to energy-driven (incident photons trapped) at this point, in 
a switch analogous to the turnoff of Gompton cooling in the wind-driven case. 
There have so far been no suggestions of how this might happen, but the physics 
of dust opacity is sufficiently complex that this is perhaps not surprising. 
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8 SMALL vs LARGE SCALE FEEDBACK 

We have shown that UFO winds are very common in AGN, despite quite restrictive 
conditions on their observability. They provide an obvious way for the central 
supermassive black hole to communicate its presence to its host. This in turn 
suggests ways of understanding both the SMBH - galaxy scaling relations, and 
the need to expel gas from the galaxy spheroid to terminate star formation. It is 
clear that this AGN feedback must operate at times on small scales, and at others 
on large scales. Our discussion of feedback points to a natural association between 
momentum—driving and small scales, and between energy-driving and large scales. 
Small-scale phenomena naturally explained by wind momentum-driving include 

1. Super-solar elemental abundances in AGN spectra. Wind momentum-driving 
automatically sweeps up and compresses the same gas many times before the 
black hole mass reaches Ma . Generations of massive stars forming out of the same 
swept-up gas can repeatedly enrich the gas close to the SMBH with nuclear- 
processed material before the Ma mass is reached, and momentum-driving changes 
to energy-driving. 

2. Dark matter cusp removal. The same repeated sweeping—up of a gas mass 
comparable to the SMBH mass, followed by fallback, has a strong tendency to 
weaken dark matter cusps. Because the baryonic mass involved is much larger, 
this is a more powerful version of the mechanism invoked by Pontzen & Governato 
(2012) (see also Garrison-Kimmel et al. 2013), who considered supernovae near 
the SMBH. 

3. Quiescence of AGN hosts. Most AGN hosts do not show dramatically ele¬ 
vated star formation in the central regions of their galaxy discs, or so far much 
evidence for high-speed (~ 1000 kms“^ and massive (~ few IOOMq yr“^) out¬ 
flows on large scales. This is compatible with wind driving by momentum but not 
energy. 

Large-scale phenomena suggesting the action of energy-driving include 

4. Metals in the circumgalactic medium. These must be made in galaxies 
and only later expelled to make the CGM. This suggests that expulsion through 
energy-driving acts only after stellar evolution has had time to enrich a significant 
fraction of the galaxy bulge gas. 

5. Mechanical luminosities of galaxy-scale molecular outflows. These are ob¬ 
served to be close to 5% of their central AGN luminosities L, just as expected for 
energy-driving, with momenta close to the predicted 20L/c. 

6. Suppression of cosmological infall. Energy-driven outflows at large radii 
probably prevent galaxies accreting indefinitely (Costa et al. 2014). 

This list seems to favor a combination of momentum and energy driving, with 
some kind of switch between them. The suppression of inverse Compton shock 
cooling at the point when the black hole mass reaches Ma appears promising, but 
requires further work on how observable the cooling is, and the possibility of two- 
fluid effects. Radiation driving on dust could produce similar behavior, although 
the physics controlling the required switch between momentum and energy driving 
is so far unexplained. 

It is worth stressing that even a detailed understanding of the dual role of 
AGN feedback in establishing both the SMBH - host scaling relations and the 
quenching of star formation would solve only half of the problem. For a full picture 
of how black holes and galaxies influence each other we need to know what physical 
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mechanism can produce a supply of gas with so little angular momentum that 
much of it can accrete on to the central supermassive black hole within a few 
Salpeter times (see Section 3.1, and equation 12011 . We saw in Section 1.2 that the 
hole’s gravity is far too weak to influence the galaxy on the mass scale needed 
for this. Only feedback can do this, perhaps suggesting that SMBH feedback may 
ultimately cause SMBH feeding (cf Dehnen & King 2013). 
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